ADA090516 


University  of  Southern  California 
Los  Angeles,  California,  90007 


.  If  C'  ■*?<> 
M*  42/ 


cu 

O 

CJ> 

<uu' 


The  views  and  conclusions  contained  in  this  document 
are  those  of  the  authors  and  should  not  be  interpreted 
as  necessarily  representing  the  official  policies, 
either  expressed  or  implied,  of  the  Defense  Advanced 
Research  Projects  Agency  or  the  US  Government. 


DISTRIBUTION  STATEMENT  A 

Approved  for  public  release; 
Distribution  Unlimited 


80  9  10  00" 


Summary 


The  purpose  of  this  project  is  the  development  of  semi-emp_r 
and  theoretical  methods  of  obtaining  reasonably  accurate  potent. .ai 
curves  for  diatomic  excimer  systems.  These  curves  are  of  importune' 
for  applications  to  fluorescence  spectra,  stimulated  emission 
coefficients,  excimer  laser  development  and  electron  impact  exci¬ 
tations.  We  have  been  mainly  interested  in  the  development  of  a  new 
effective  potential  method  which  obtains  potential  curves  by 
calculating  the  energy  levels  of  one  atom  in  the  external  field  of 
the  other  atom.  In  order  to  be  able  to  compare  our  new  method  to 
other  well  established  methods,  we  have  first  done  a  Cl  calculation 
for  GaKr  and  applied  the  Gordon-Kim  method  to  GaKr.  This  provided 
some  knowledge  of  the  effect  of  a  Kr  atom  on  another  atom  which 
might  be  used  later  as  input  into  the  effective  potential  method. 

The  basic  idea  of  the  new  method  is  as  follows:  In  all  low 
lying  states  of  a  metal-rare  gas  system  one  atom  (the  rare  gas  atom 
which  has  a  high  first  excitation  energy)  is  asymptotically  in  its 
ground  state  and  the  various  molecular  states  correspond  to  the 
various  low  lying  excited  states  of  the  metal  atom.  Therefore  the 
basic  assumption  of  our  method  is,  that  it  is  possible  to  calculate 
the  molecular  energy  levels  as  the  energy  levels  of  a  single  metal 
atom  in  an  external  field  which  represents  the  rare  gase  atom.  This 
external  field  which  represents  the  second  atom  is  our  effective 
potential  The  general  form  of  this  potential  can  be  derived  from 
the  many  body  Greens  function.  Phenomenological  treatments  of  the 
many  body  theory  leads  to  a  form  with  free  parameters.  Our  idea  is 
to  determine  these  parameters  by  fitting  one  molecular  potential 
curve,  calculated  by  our  new  method,  to  a  known  curve.  We  hope  that 


the  effective  potential  resulting  from  one  particular  state  describes 

the  effect  of  the  rare  gas  atom  also  in  other  molecular  st  _es  and 

even  in  other  molecules,  i.e.  independently  of  the  type  and  stat-. 

of  the  other  atom.  The  big  problem  in  this  new  method  is,  that  new 

types  of  integrals  have  occured,  which  have  not  appeared  in  quantum 

chemistry  before  and  for  which  no  standard  computer  programs  exist. 

Therefore,  we  spent  most  of  the  time  developing  and  testing  new 

computer  programs.  By  now  the  programs  seem  to  work  properly  but 

they  are  still  very  slow  and  at  the  moment  we  are  still  testing 

the  whole  method  on  a  simple  system,  on  LiHe.  We  have  tried  to 

parametrize  the  effective  potential  of  the  He  ground  state  by 

2 

fitting  the  resulting  molecular  potential  for  the  rr  state  of  HeLi 
to  the  known  potential  curve  of  this  state.  Probably  we  did  not 
use  the  best  fitting  strategy  so  far  and  therefore  we  did  not  get 
a  satisfactory  fit  of  the  molecular  potential.  It  seems  to  be 
necessary  to  repeat  the  fitting  procedure  with  a  more  efficient 
strategy,  which  we  consider  in  the  text.  The  next  step  in  the 

program,  after  the  fit  of  the  molecular  potential  curve  for  one 
state,  is  to  calculate  the  molecular  potential  curves  for  other 
states  of  HeLi  using  the  same  effective  potential  to  represent  the 
He  atom.  The  result  of  this  test  will  show  if  the  whole  idea  works 
or  not.  We  are  continuing  work  on  this  idea  even  without  the 
support  of  ARPA-ONR  as  we  believe  it  to  be  a  valuable  new  method. 
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1 . 1  Introduction 

The  general  purpose  of  this  project  is  to  continue  the  developn-' 

and  implementation  of  semi-empirical  and  theoretical  methods  of 

obtaining  reasonably  accurate  potential  curves  for  diatomic  excimer 

systems.1  Once  the  potential  curves  have  been  obtained,  they  can 

2 

be  used  to  predict  fluorescence  spectra  ,  and  to  compute  laser  gain, 
stimulated  emission  coefficients3  and  total  cross  sections  for 
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excitation. by  electron  impact.  While  highly  accurate  potential 
curves  are  always  a  goal,  our  immediate  goal  is  the  relatively  rapid 
production  of  potential  curves  that  are  accurate  enough  to  guide 
experimentalists  in  choosing  or  rejecting  possible  laser  systems 
on  the  basis  of  inexpensive  theoretical  calculations  rather  than 
on  the  basis  of  expensive  and  time-consuming  experiments. 

We  are  particularly  interested  in  developing  methods  that  in 
the  future  would  be  applicable  to  excimer  systems  because  of  the 
current  emphasis  on  these  systems  as  candidates  for  efficients, 
high-power  visible  and  ultraviolet  lasers.1  The  specific  systems 
that  motivated  this  proposal  were  the  group  IIA-,IIB-  and  IIIA- 
rare  gas  systems  (such  as  MgXe  and  TlXe) ,  alkali-group  IIB  systems 
and  Thallium-group  IIB  systems.  Our  immediate  goal  is  to  develop 
a  new  efficient  method  to  calculate  potential  curves  for  non-bounding 
diatomic  systems  in  which  the  asymptotic  atoms  have  greatly  different 
ionization  potentials  (as  do  most  metal-rare  gas  systems).  We  will 
then  view  the  lowest  several  states  of  the  molecule  as  roughly 
represented  by  a  specific  state  of  the  metal  in  the  effective 
"external"  field  of  the  high  ionization  potential  (i.e.  here  a 
rare  gas)  atom. 


The  next  step  will  be  to  parametrize  this  potential  (and  the 
effective  electron  potentials  that  arise  therein)  using  a  form 
suggested  by  field  theory  and  data  from  any  known  molecular  state 
of  any  metal-rare  gas  system  that  goes  at  large  internuclear 
distance  to  the  ground  state  of  the  rare  gas.  This  effective 
potential  will  then  be  used  to  describe  unknown  states  of  the  same 
or  different  metals  by  viewing  their  calculation  as  one  of  the 
metal  atom  in  the  field  of  this  potential.  If  this  aim  is  attained 
a  basis  empirical  transfer  of  information  from  known  to  unknown 
systems  and  states  will  be  achieved  with  an  in  principle  accompaning 
simplification  in  calculations.  This  in  turn  would  allow  us  to  go 
to  bigger  systems. 

Although  extensive  work  has  been  done  on  calculating  potential 
curves,  the  currently  available  methods  are  inadequate  for  large 
systems  such  as  TlXe.  While  self-consistent  field  (SCF)  and 
configuration  interaction  (Cl)  programs  are  the  most  accurate,  at 
this  time  those  computer  codes  cannot  handle  molecular  systems 
with  f  electrons,  such  as  TlXe. 

Now,  simplified  methods,  that  are  alternatives  to  Cl  do  exist: 

The  Gordon-Kim^  (GK)  method  is  available  but  has  not  been  adequately 
tested  for  large  systems  with  open  shells.  We  believe  that  our 
method,  theoretically  is  better  based  than  this  method,  as  the 
reader  can  judge  for  himself  in  later  sections,  and  it  takes  in 
such  effects  as  mutual  atomic  polarization  that  the  GK  leaves  out. 

The  most  promising  method  for  treating  big  systems  rigorously 

is  the  pseudopotential  method.  The  pseudopotential  method**  is 
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currently  being  extended  to  large  systems.  ' 
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The  new  method  we  are  developing,  the  effective  potential 

method,  is  essentially  a  means  of  calculating  the  potential  curv 

for  a  diatomic  system  AB  (where  A  is  a  closed  shell  atom)  by  doing 

an  atomic  SCF  or  Cl  calculation  on  B  in  which  A  is  treated  as  an 

external  (effective)  potential.  The  effective  potential  is  the  sun- 

of  the  Hartree-Fock  potential  for  A  and  a  polarization  potential . 

The  polarization  potential  consists  of  a  one-particle  part  describing 

the  polarization  of  A  by  an  external  charged  particle  and  a  two- 

particle  part  which  corrects  for  the  "depolarization"  due  to  the 

presence  of  the  other  electrons  and  nucleus  of  3.  The  parameter  of 

the  effective  potential  will  be  obtained  by  using  perturbation 

theory  and  a  known  potential  curve  of  AX  (X  is  any  atom).  If  X  is 

not  highly  polarized,  a  first  order  approximation  to  the  effective 
A 

potential  U  may  be  obtained  by  fitting  the  parameters  to 

VX,e(R)  =  ^dr  UA(R,i)px(0 

A 

where  p  is  the  density  of  atom  X  and  Vv,e  is  the  electronic  part 

A  A 

of  the  known  potential  curve  for  AX.  The  theory  of  the  effective 
potential  and  the  perturbation  treatment  of  it  are  discussed  in 
more  detail  in  §4. 

This  effective  potential  method  has  several  advantages  over  the  Gordon-KIm 
and  pseudopotentlal  methods.  Because  the  effective  potential  allows  for  polar¬ 
ization  of  A  and  the  SCF  calculation  allows  for  polarization  of  B,  this  method  should 
be  more  accurate  than  GK.  At  the  same  time,  LS  coupling  and  relativistic  effects  can 
Included  as  easily  In  our  method:  the  treatment  of  LS  coupling  would  be  the  same  and 
relativistic  effects  can  be  Included  In  the  core  potentials.  In  comparing  our 
effective  potential  method  to  pseudopotential  methods,  the  following  points  should  be 
noted.  Methods  which  use  only  a  pseudopotential  to  replace  the  valence-core  Inter- 
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actions  do  not  include  polarization  of  the  core.  In  contrast,  our  methcc  'dudes 
core  polarization;  consequently,  the  core-core  interaction  is  correct  ?t  i 
distances  so  the  correct  van  der  Waals  interaction  will  be  obtained.  Becau"?  w-~ 
take  advantage  of  the  fact  that  there  is  no  real  bonding  in  the  systems  we  are 
studying,  our  calculations  will  be  essentially  atomic  calculations  and  will  include 
fewer  electrons  explicitly  than  a  pseudopotential  calculation.  Use  of  a  sem*- 
empirical  potential  which  is  parameterized  to  molecular  data  should  result  in  better 
molecular  potential  curves  than  use  of  semi -empirical  potentials  parameterized  to 
atomic  data.  In  addition,  use  of  molecular  data  will  allow  us  to  obtain  effective 
potentials  for  electrons  outside  of  rare  gas  atoms,  and  for  electrons  outside  of 
highly  charged  cores;  these  are  two  cases  in  which  parameterization  to  atomic 
data  would  be  extremely  difficult. 

Since  the  effective  potential  method  is  semi -empirical  and  since  the  potential 
parameters  will  be  fit  to  molecular  data,  some  potential  curves  are  needed.  Basically, 
there  are  two  options  available  for  parameterizing  the  effective  potential  of  A  in 
order  to  calculate  potential  curves  for  AB:  1)  The  effective  potential  for  A. can 
be  fit  to  known  potential  curves  for  XA  (where  X  is  any  atom)  or  2)  The  effective 
potential  for  A  can  be  fit  to  whatever  ground  and  excited  state  curves  for  AB  are 
known.  The  second  option  recognizes  the  facts  that  it  is  easier  to  obtain  ground 
state  curves  than  excited  state  curves  by  other  methods  and  that  both  are  needed  to 
pick  possible  laser  systems.  The  curves  needed  for  fitting  can  be  obtained  by  a 
variety  of  techniques  varying  from  rigorous  Cl  calculations  to  analysis  of  experimental 
data.  As  examples  of  available  techniques,  consider  the  following:  When  experimental 
band  spectra  data  exist,  the  quasi-static  theory  of  line  shapes  could  be  used  to 

g 

extract  potentials.  Atomic  beam  data  (such  as  that  supplied  by  Yuan  Lee  In 

Berkeley)  could  also  be  used  to  fit  potential  forms.10  The  RKR  method11  could  be 

used  when,  as  is  rarely  the  case  in  excimers,  vibrational  data  exists.  For  positive- 

12 

negative  Ion  systems  Rlttner  potentials  can  be  fit  to  available  data.  For  systems 
with  closed-shell  separated  atoms  and  ions  or  for  certain  types  of  open-shell 
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systems,  the  Gordon-Kim  method"*  can  be  used  to  predict  short  range- 
potentials,  and  long  range  corrections  can  be  added  to  this  model. 

Short  range  repulsive  potentials  for  closed-shell  atoms  and  ions 
can  also  be  fitted,  correlated  and  interpolated  using:  ( a t  the 

*  delta-function  model  for  homonuclear  systems10;  (b)  the  delta- 
function  and  distortion  model  with  combination  rules  for  hetero- 

V  nuclear  systems10;  (c)  the  related  softness  parameter  methods  that 

1  3 

require  atomic  Hartree-Fock  charge  densities. 

Now  after  listing  the  advantages  of  our  method  we  must  admit 
to  slow  progress  and  a  lack  of  success  in  so  far  completing  our 
objectives  due  to  two  problems.  The  first  is  that  our  effective 
potential  leads  to  new  two  center  (atom-external  potential)  integrals 
such  as  those  representing  the  interaction  with  the  part  of  the 
effective  potential  representing  the  polarized  rare  gas  (see  luLer 
sections  for  details) .  We  have  developed  methods  and  programs  feu 
treating  these  integrals  but  unfortunately  the  programs  are  still 
very  consuming  of  computer  time.  These  integrals  have  never  appeared 
in  quantum  chemistry  before  and  we  had  to  develop  new  methods  of 
of  doing  them.  We  do  not  believe  our  methods  are  presently  optimum. 

Big  savings  in  computer  time  can  probably  be  achieved  with  a  few 

*  years  more  efforts  ("years"  may  seem  long  but  not  when  one  considers 
the  thirty  or  so  years  it  took  to  reach  the  state  of  the  art  in 

,  present  molecular  programs).  Our  strategic  problem  was  that  it  was 

senseless  to  develop  speedy  packages  to  use  in  a  method  that  might 
not  work.  Hence  we  used  our  slow  packages  and  went  on  to  the  second 
problem,  that  of  parametrizing  the  effective  potential.  In  a  sense 
the  problems  are  linked,  in  that  the  trial  and  error  fittings 

»  i' 

necessary  to  parametrize  are  severly  limited  in  application  by  the 
computer  cost  for  doing  integrals  as  parameters  are  varied.  The 
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second  problem  is  therefore  a  so  far  less  than  ideal  but  per  ■'ps 
not  unacceptable  fit  of  the  parameters  in  the  effective  ootential 
to  the  known  curves.  We  will  return  to  this  problem  in  section  4.V. 

Now  of  the  outset  of  our  work  we  decided  that  we  would  Like 
as  input  from  theory  some  accurate  GaKr  curves  so  that  we  could 
get  on  effective  potential  for  Kr  to  use  in  InKr  or  TlKr  computations. 
Hence  we  split  our  work  in  three  parts  -  the  first  part  reviewed  in 
§2  was  to  do  accurate  Cl  calculations  on  GaKr. 

§3  is  an  application  of  the  Gordon-Kim  method  to  GaKr  and  a 
comparison  with  the  results  from  §2. 

The  third  part  was  to  develop  and  test  our  theory  in  §4.  As  it 
turned  out  the  third  part  is  still  not  able  to  take  advantage  of 
the  first  part  since  even  for  Li-He  we  have  not,  because  of  the 
above  two  mentioned  problems,  achieved  success.  In  any  case  we  will 
report  here  our  progress  ommitting  the  large  sections  of  efforts 
that  lead  to  even  inferior  methods  of  integration  and  parametri/ation 
and  reporting  only  the  most  promising  methods.  At  the  writing  of 
this  report  we  are  close  to  a  critical  test  calculation  of  Lille, 
having  achieved  what  may  be  a  satisfactory  (how  satisfactory  it  is 
future  results  will  show)  calculation  of  the  effective  Helium 
potential.  Time  and  personnel  changes  have  prevented  us  from  making 
this  last  step  which  we  shall  certainly  do  in  the  near  future  even 
without  further  support.  If  we  are  successful  we  shall  return  to 
the  problems  of  speeding  up  our  integral  computation  methods  and 
routines.  This  will  prepare  the  method  to  move  on  to  bigger  systems 
and  get  us  to  where  we  had  hoped  to  be  now,  at  the  end  of  our  contract 
period.  Of  course  if  the  calculations  are  not  successful  theoretical 
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reasons  need  bo  sought  and  the  promise  of  the  method  would  bo  in 
question.  While  doing  this  work  Dr.  Taylor  and  Dr.  Jung  had  a 
discussion  about  electron  scattering  in  a  strong  laser  field 
and  wrote  a  paper  about  this  subject .  A  copy  is  included  in 
Appendix  E.  This  work  was  not  in  the  original  grant  but  as  many 
things  in  science  it  just  happened  in  discussion.  No  computing 
fonds  from  the  grant  were  used  for  this  project. 
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§2 

Flectronic  States  of  GaKr  :  Ab  initio  calculations  of  a  prototype  for 
TlKr 


2.1  Introduction 

Among  the  metal-rare  gas  eximers,  the  TIXe  system  is  believed  to 

be  an  excellent  candidate  for  an  efficient  visible,  high  power,  tunable 

laser.1  Since  <^b  initio  calculations  on  this  system  are  beyond  the  scope 

of  present  computer  programs,  we  present  here  a  configuration  interaction 

■Cl )  calculation  on  GaKr  ,  which  is  the  largest  group  IIIB-rare  gas  system 

for  which  Cl  calculations  can  be  done.  We  use  the  calculated  GaKr  curves 

to  model  the  potential  curves  for  InKr  and  TlKr  .  Although  our  model 

does  not  allow  further  extrapolation  from  TlKr  to  TIXe  ,  it  is  hoped 

that  these  calculations  will  yield  some  insight  into  the  properties  of  the 

2 

TIXe  eximer.  In  addition,  Gallagher  has  recently  raised  the  possibility 
of  using  GaXe  as  a  laser  if  the  Ga  can  be  obtained  from  dissociation 
of  Gal^  •  Consequently  the  GaKr  curves  should  also  be  of  intrinsic 
i nterest. 


In  this  paper,  the  Cl  calculations  on  GaKr  are  presented  alonq  with 
the  model  calculations  on  InKr  and  TlKr  .  These  potential  curves  are 
used  as  the  basis  of  a  classical  calculation  of  the  emission  and  absorption 
coefficients  for  these  systems. 
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2.  II  Details  of  the  calculation 

A.  Basis  set 

The  calculations  use  (14sllp6d)  primitive  Gaussian  bases  for  gallium 

3 

and  krypton  as  a  starting  point.  The  core  orbitals  (Is  ,2s ,2p,3s  ,3p)  are 
singly  contracted  to  the  Hartree-Fock  atomic  orbitals  while  the  valence 
orbitals  (4s  and  4p)  and  the  3d  orbitals  are  each  described  by  two  contracted 
functions (see  Table  I).  The  resulting  (5s4p2d/5s4p2d)  contracted  bases 

4 

are  constructed  using  the  general  contraction  scheme  of  Raffenetti . 

These  basis  sets  are  extended  to  include  polarisation  by  adding  two 
diffuse  s  functions  (5s  and  5s')  and  a  diffuse  p  function  (5p)  to 
describe  the  lowest  Rydburg  orbitals.  The  exponents  for  these  orbitals 
(c(5s)  =  .026,  libs')  -  .011,  t(5p)  =.01)  are  obtained  from  atomic  calcu¬ 
lations  on  the  excited  states  of  Ga  .  The  final  basis  set  thus  consists 
of  a  ( 16s 1 2p6d )  primitive  basis  contracted  to  [7s5p2d]  for  Ga  and  a 
(14sl lp6d)  primitive  basis  contracted  to  [5s4p2d]  for  Kr  . 

« 

;  B.  SCF  calculation 

! 

|  ! 

The  starting  point  for  the  Cl  calculation  is  a  Hartree-Fock  calculation 
2  + 

on  the  t.  state 

I 

I  ....  13n214nZ7n415o216n 


-  11  - 

The  inner  core  molecular orbi tals  (twelve  n  ,  twelve  n  and  four  4) 
are  held  doubly  occupied  from  this  point  on  and  are  replaced  by  the 
rigorous  nonlocal  Hartree-Fock  potential. 

*core-,N*1E  <ZJI  -  Ki> 
i=core 

With  the  core  orbitals  removed  from  consideration,  it  is  convenient 
to  renumber  the  valence  orbitals  so  that  the  Hartree-Fock  configuration 
is  written  as 

,  2,  2-  2.  ,  4 

io  2k1  3a 

At  large  R  the  correspondence  for  the  valence  orbitals  is 

la  *  4s Kr 

2n  »  4pKr 

3a  -*■  4s  Ga 

4o  *•  4pGa 

In  -►  4pKr 
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In  addition  to  the  valence  orbitals,  nine  a  ,  six  *  and  four  6 
virtual  orbitals  are  used  in  the  Cl  calculations.  The  lowest  virtual  orbitals 
(5o ,  6a,  7a,  2v ,  3n),  which  correspond  to  the  Ga  5s,  5s‘,  4p  and  5p 
atomic  orbitals  for  large  internuclear  separations,  are  obtained  by  the 
improved  virtual  orbital  (IVO)  procedure.^  The  IVO  orbitals  are  obtained 
by  removing  the  electron  from  the  4a  (valence)  orbital  of  the  above 
configuration  and  calculating  the  virtual  orbitals  for  the  (N-l )-electron 
Hami 1  torn  an. 


C.  Cl  calculations 

Full  polarization  Cl  (POL-CI)  calculations^  which  provide  a  balanced 

description  of  all  states  of  the  4s-4p  manifold  were  carried  out.  A  set 

of  reference  configurations  was  chosen  (see  Table  II)  to  describe  the  dominant 

configurations  for  the  2r.+  and  2n  states  of  the  molecule  and  the 
1  + 

state  of  the  ion. 

The  full  POL-CI  calculations  include  all  (1+2)  electron  excitations 
relative  to  each  reference  configuration  subject  to  the  restrictions  that 
no  more  than  one  electron  occupy  the  Rydberg  5 a  orbital  and  no  more  than 
one  electron  occupy  any  virtual  orbital  (6o,  2n  etc.).  This  results  in 

P  x 

764  spatial  and  2314  spin  configurations  for  the  T.  states,  556  spatial 
and  1565  spin  configurations  for  the  n  states  and  368  spatial  and  558 
spin  configurations  for  the  ^  state.  A  total  of  15a  ,  14*  and  26 
occupied  and  virtual  orbitals  are  used  in  the  POL-CI  calculations. 
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JRECHttwa  pack  bumc-nct  nuts 


The  potential  energy  curves  (Fig.  2  and  Table  III)  from  the  POL-CI 

calculations  generally  follow  the  behavior  predicted  from  these  theoretical 

considerations.  These  calculations  are  not,  however,  designed  to  treat 

long-range  dispersion  forces.  Many  excitations  which  contribute  to  a 

C^r  ^  attraction  are  not  included  in  the  wavefunction.  For  these  reasons, 

one  would  expect  deeper  wells  in  all  the  potential  curves  including  the 

2 

two  lowest,  generally  repulsive,  states.  As  predicted, the  1  n  state 

2  +  ?  + 
is  less  repulsive  than  the  1  i  state.  The  bound  excited  state  ?  ' 

has  a  minima  at  6.36  a.u.  about  3.0  eV  above  the  ground  state  curve,  while 

the  minima  in  the  ion  curve  is  at  6.28  a.u.  and  5.6  eV  above  the  ground 

state.  The  well  depths  are  compared  with  those  obtained  by  Gallagher11 

in  Table  IV,  and  as  expected,  Gall^her's  wells  are  deeper. 

The  dipole  moments  of  these  states  and  the  transition  moments  between 

2  + 

T.  and  the  lower  states  have  also  been  calculated  and  are  given  in 
Table  V  and  Figs.  3  and  4. 
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B.  Electronic  states  including  spin-orbit  coupling 

A  complete  treatment  of  the  electronic  states  of  GaKr  must  include 
the  effects  of  spin-orbit  coupling.  The  states  considered  here,  which 
dissociate  to  the  closed-shell  ground  state  of  Kr  and  an  open-shell  state 
of  Ga  or  Ga*  are  influenced  only  by  the  soin-orbit  matrix  elements  oi 
the  open- she  11  atom. 

Following  the  procedure  used  previously, we  have  adopted  a 
simple  model  for  including  the  effects  of  spin-orbit  coupling  on  the 
calculated  potential  energy  curves  and  wavefunctions .  The  experimental 
spin-orbit  parameters  for  the  open-shell  atom  (Ga  and  later  In  ,  Tl)  are 
used  to  determine  the  matrix  elements  of  the  spin-orbit  interaction,  H  , 
coup  I i ng  the  molecular  states  at  infinite  separation  and  these  matrix 
elements  are  assumed  to  be  independent  of  internuclear  distance.  The 
resulting  spin-orbit  matrix  Hsq  is  added  to  the  diagonal  matrix  of 
electronic  energies  He|ec(^)  '  <5^jE^(R)  : 

H<"> '  W">  *  % 

and  the  total  matrix  H  Is  then  diagonalized  at  each  internuclear  distance 
R  .  Thus,  in  addition  to  the  assumption  that  the  spin-orbit  matrix  elements 
do  not  change  as  functions  of  R  ,  this  model  assumes  that  only  one-center 
terms  need  be  included  so  that  only  the  spin-orbit  coupling  on  Ga  is 
important.  He  would  expect  this  procedure  to  provide  reasonable  results 
to  the  extent  that  the  molecular  states  retain  the  identity  of  the  atomic 
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states  from  which  they  are  formed.  The  spin-orbit  matrices  are  given 

in  Table  VI  along  with  the  atomic  parameters  used  in  these  calculations. 

2  2 

The  parameter  x  is  chosen  so  that  the  atomic  an<*  Pl/2  states 

have  energies  of  +  x  and  -2x  respectively. 

We  shall  label  the  molecular  states  using  the  convention  of  Hund’s 

case  (c)  where  i.  ,  the  projection  of  total  angular  momentum  along  the 

molecular  axis,  is  the  only  good  quantum  number,  si  is  defined  as 

-  ’  ♦  Sj  ,  where  i  and  S^  are  the  orbital  and  spin  angular  momentum 

2  + 

projection,  respectively.  The  molecular  r.  states  have  only  a 

2 

.  -  1/2  component,  while  the  r  states  yield  a  sj  -  3/2  and  ft  =  1/2 
state  the  states  in  the  ft  representation  are  labeled  according  to 
increasing  energy  by  a  Roman  numeral.  So  the  ft  -  1/2  states  are  designated 
as  M/2,111/2  and  the  Si  -  3/2  states  are  13/2,  II  3/2  . 

The  coupled  states  are  expressed  as  follows: 

111/2-  -  c,.|lV>  +  Cn|l?n> 

-II  1/2  >  =  -cn|lV>  +  Cj.|l2n> 

1 1  3/2  >  -  1 1 2 n » 

The  spin-orbit  coefficient  CE  and  Cn  are  given  In  Table  VII.  The 
potential  curves  and  transition  moments  for  GaKr  Including  spin-orbit 
coupling  are  given  In  Table  VIII  and  Figs.  5  and  6. 

V 

i 

/ 
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2. IV  Extrapolations  to  InKr  and  TlKr 

The  potential  curves  for  InKr  and  TlKr  are  modeled  on  the 

GaKr  curves.  The  lowest  excitation  energies  and  the  ionization  potentials 

for  the  series  Ga  ,  In  ,  T1  are  given  in  Table  IX.  As  can  be  seen,  this 

series  does  not  form  a  steady  progression.  In  has  a  lower  ionization 

potential  and  lower  excitation  energies  than  Ga  ,  as  expected  for  a 

heavier  atom.  However,  T1  has  a  higher  Ionization  potential  and  higher 

excitation  energies.  This  is  due  in  part  to  the  presence  of  a  filled 

12 

4f  shell  in  T1  and  the  larger  spin-orbit,  effects.  These  effects 
should  be  considered  when  extrapolating  the  GaKr  curves  to  InKr  and 
TlKr. 


To  simulate  InKr  and  TlKr  ,  the  experimental  spin-orbit  parameters 
for  In  and  T1  are  used  to  couple  the  GaKr  curves.  The  curves  are 
also  shifted  to  give  the  correct  atomic  excitation  energies  at  R  =  ~ 

(see  Table  IX).  This  procedure  should  give  at  best  a  qualitative  description 
of  states  of  InKr  and  TlKr  ,  since  the  non-spin-orbit  coupled  states 
are  expected  to  have  quantitatively  different  well  depths  and  equilibrium 
separations . 

The  effect  of  the  increasing  spin-orbit  perturbation  in  going  from 
Ga  to  In  and  T1  is  evident  in  the  calculated  curves  which  are  given 
in  Figs.  7  and  8.  Only  the  well  depths  and  positions  for  the  I  1/2  and 
II  1/2  states  are  affected  by  the  spin-orbit  coupling.  The  other  states 
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are  the  same  as  those  for  GaKr  except  that  they  have  the  correct 
asymptotic  spacing. 

The  mixing  parameters  from  the  spin-orbit  coupling  calculations  for 
T1 Kr  are  also  used  as  the  coefficients  of  the  GaKr  wave  function  to 
estimate  the  transition  moments  for  this  system.  The  TIKr  transition 
moments  are  given  in  Table  X  and  Fig.  9. 
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2.V  Absorption  and  stimulated  emission  coefficients  for  possible 

laser  transitions 

The  Interest  in  the  group  IIIB-rare  gas  systems  arises  from  the  possibility 
of  their  use  as  visible  laser  systems.  In  order  to  judge  their  usefulness  as 
lasers  it  is  convenient  to  calculate  the  absorption  k^(t)  and  stimulated 
emission  g  ( T )  coefficients.  Obtaining  quantum-mechanical  results  for  these 
quantities  would  require  a  complex  calculation  which  would  be  inconsistent  with 
the  extrapolations  used  to  obtain  the  InKr  and  TIKr  curves.  Consequently, 
we  have  used  Gallagher's  analysis  [17].  which  is  based  on  the  classical 
f rank- Condon  principle. 

In  order  to  obtain  g  and  k  ,  the  Cl  curves  are  first  fit  by  Morse 
potentials.  The  parameters  for  these  Morse  potentials  are  given  in  Table  XI. 
These  parameters  can  then  be  used  In  Gallagher's  equations,  along  with  the 
atomic  transition  rate,  to  obtain  absorption  and  stimulated  emission  coefficients 
for  pressure  and  excitation  conditions  of  interest  to  experimentalists.  We 
have  calculated  these  coefficients  for  two  different  types  of  conditions. 

The  high  temperature  results  correspond  to  the  case  where  the  concentration 
of  the  metal  is  obtained  from  the  vapor  pressure  of  the  metal  itself,  while 
the  low  temperature  results  correspond  to  obtaining  the  required  concentration 
of  the  metal  from  vaporization  of  MI3  (M  =  Ga,  In,  or  Tl).  This  latter 
condition  has  been  suggested  by  Gallagher  as  a  possible  means  of  obtaining 
high  concentrations  of  the  metal  at  low  temperatures.  In  both  cases  the 

densities  used  are  1020/cm3  for  Kr  ,  1016/cm3  =  3fM  2P,  =  l  5[m  2P  ] 

14  3  2  '  C  '  3/2 

and  2  x  10  /cmJ  *  [M  •  The  resulting  absorption  and  stimulated 
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Missfon  coefficients  for  GaKr  ,  InKr  end  TIKr  are  given  In  Figures  10 
to  15. 

A  paper  has  been  published  about  the  work  presented  in  §2. 

See  Appendix  A. 
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Table  I.  Gaussian  exponents  and  contraction  coefficients 


Is 

2s 

457600. 

.  000222 

-.000069 

68470. 

.001732 

-.000535 

15590. 

.008952 

-.002814 

4450. 

.035874 

-.011275 

1472. 

.114000 

-.038495 

541.3 

.274138 

-.100714 

214.8 

.414793 

-.211832 

88.81 

.275395 

-.175448 

27.18 

.029561 

.479840 

11.54 

-.006815 

.634145 

3.303 

. 002253 

.069592 

1.334 

-.001017 

-.012299 

.1947 

.000251 

.002774 

.07158 

0.0 

no 

3s 

4s 

4s' 

. 000026 

-.000006 

0.0 

. 000205 

-.000048 

0.0 

.001070 

-.00024/ 

0.0 

.004337 

-.001007 

0.0 

.014707 

-.003399 

0.0 

.039748 

-.009279 

0.0 

. 084475 

-.019587 

0.0 

.079654 

-.019104 

0.0 

-.291821 

.072753 

0.0 

-.527118 

.134137 

0.0 

.583707 

-.181778 

0.0 

.674103 

-.358241 

0.0 

.028077 

.615164 

0.0 

0.0 

0.0 

1.0 

5s 

■026  l.o 

•Oil  0.0 


5s' 

0.0 

1.0 


3274. 

765.4 

241.6 

89.39 

36.36 

15.60 

6.472 

2.748 

1.090 

.2202 

.06130 

.01 


2p 

.001513 
.013070 
.067263 
.219542 
.421107 
.376515 
.089425 
-.000502 
.001 761 
-.000247 
0.0 

0.0 


3p 

-.000576 
-.004981 
-.026421 
-.089529 
-.186734 
-.144494 
.258956 
.570187 
. 325305 
.016563 
0.0 

0.0 


4p  4p’  5p 


. 000094 
.000800 
.004337 
.014443 
.031377 
.021501 
-.046233 
-.125293 
-.045636 
.452811 
0.0 

0.0 


0.0  0.0 

0-0  0.0 

0-0  O.o 

0.0  0.0 

0.0  o.o 

0.0  o.o 

0-0  o.o 

0.0  0.0 

0.0  0.0 

0.0  0.0 

1-0  o.o 

0.0  1.0 


3d 

59.66  .031949 

17.10  .163546 

6.030  .367457 

2.171  .456851 

.6844  .305161 

.160  0.0 


4d 

0.0 

0.0 

0.0 

0.0 

0.0 

1.0 
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Table  I.  Gaussian  exponents  and  contraction  coefficients 


Exponents 

Contraction 

coefficients 

Krypton  atom 

Is 

2s 

3s 

4s 

4s' 

605700. 

.000231 

-.000073 

.000029 

- . 000009 

0.0 

90300. 

.001755 

-.000551 

.000221 

-.000070 

0.0 

20920. 

.009076 

-.002894 

.001159 

-.000369 

0.0 

5889. 

.036990 

-.011834 

.004781 

- . 001 522 

0.0 

1950. 

.116154 

-.039826 

.016056 

-.005118 

0.0 

718.2 

.278401 

-.104801 

.043454 

-.013886 

0.0 

285.4 

.415746 

-.217093 

.091899 

-.029537 

0.0 

118.6 

.267204 

-.175562 

.083789 

-.027309 

0.0 

38. 16 

.027870 

.471395 

-.303023 

.103498 

0.0 

16.45 

-.005998 

.636794 

-.570620 

.208810 

0.0 

5.211 

.002217 

.082255 

.501751 

-.235737 

0.0 

2.291 

-.001092 

-.014138 

.760483 

-.553570 

0.0 

.4837 

.000306 

.003289 

.044857 

.701123 

0.0 

.  1855 

0.0 

0.0 

0.0 

0.0 

1.0 

2p 

3p 

4p 

4p ' 

4678. 

.001392 

-.000569 

.000156 

0.0 

1120. 

.011666 

-.004777 

.001286 

0.0 

357.  1 

. 060858 

-.025631 

.007059 

0.0 

131.4 

.210040 

-.092159 

.024990 

0.0 

52.86 

.421000 

-.200936 

.056870 

0.0 

22.  70 

.383515 

-.160784 

.040225 

0.0 

9.547 

.097383 

.267789 

-.084756 

0.0 

4.  167 

-.001087 

.585908 

-.240291 

0.0 

1.811 

.002209 

.291397 

-.038636 

0.0 

.5337 

-.000509 

.015484 

.599154 

0.0 

.1654 

0.0 

0.0 

0.0 

1.0 

3d 

4d 

125.6 

.019168 

0.0 

33.31 

. 125638 

0.0 

12.15 

.366069 

0.0 

4.350 

.502482 

0.0 

1 .494 

.264377 

0.0 

.35 

0.0 

1.0 
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Table 


.  Reference  Configurations 

2  + 

l  states 

,  ,  2,  2,  2.  .  4 

1  la  2a  3a  4g1tt 

2  1 a22a2 3o25al n4 

2  2  2  4 

3  la  2a  3a  6a In 

2n  states 

1  lCT22o23n2U427T 

2  la22a23a2l7T43n 

state  (GaKr+) 
1  la22a23a21ir4 
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Table  III. 

POL-CI  calculations  on  the  low-lying 
and  the  ground  state  of  GaKr+  .  All 
relative  to  -4674.  hartrees. 

states  of  GaKr 
energies  are 

R 

iV 

2V 

12J1 

iV 

do 

-1.200042 

-1.095050 

-1.198917 

-0.992860 

15.00 

-1.200121 

-1.095025 

-1.198978 

-0.993044 

10.00 

-1.200290 

-1.094061 

-1.199328 

-0.994018 

8.00 

-1.198650 

-1.094202 

-1 .200204 

-0.996703 

7.00 

-1.193963 

-1.095454 

-1.200450 

-0.999505 

6.00 

-1.179381 

-1 .096768 

-1.197718 

-1.002018 

5.00 

-1.135164 

-1 .088552 

-1.179245 

-0.992868 

4.50 

-1.086883 

-1 .064404 

-1.148099 

-0.967583 

4.00 

-1.018697 

-0.992381 

-1.079327 

-0.906827 

3.75 

-0.968656 

-0.933786 

-1.025394 

-0.856247 
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Table  IV.  Potential  well  depths 


Mole 

cule 

State 

Cl 

R 

AE(eV) 

Morse 

R 

Fit 

AE(eV) 

Gallagher 

R  AE(eV) 

GaKr 

Il/2(Xl/2) 

^.018 

7.55 

.021 

r3/2(X3/2^ 

^-.040 

7.18 

.040 

Hl/2(Al/2* 

MO 

.006 

10.15 

.00642 

iii1/2(b2e1/2) 

6.36 

.064 

6.26 

.080 

GaKr+  IQ 

6.28 

.26 

5.95 

.252 

InKr 

1 1  /2 

•vfl 

^.013 

7.86 

.00778 

I3/2a 

1 1 1  /  2 

MO 

^.006 

9.98 

.0064 

III  a 

1  M/2 

'oa 

TIKr 

1 1  /  2 

^8 

M010 

8.5 

.012 

7.01 

.024 

1 3/2a 

6.58 

.062 

ni/2 

MO 

%.006 

9.88 

.064 

ini/2a 
,  a 

6.09 

.107 

asame  as  GaKr 


27 


Table  V. 

Dipole  and 

transition  moments  for  the  low-lying  states  of  GaKr 

R 

2 

XI 

12E+  22I+  12E+-X2n  22Z+-X2n  22E+-12E 

15.00 

0.00598 

0.00156 

0.07558 

0.00183 

-1.29167 

-1.31011 

10.00 

0.02605 

0.02048 

0.56872 

0.01631 

-1.28976 

-1.29762 

8.00 

0.09152 

0.14168 

0.89914 

0.04761 

-1.27767 

-1 .23098 

7.00 

0.18861 

0.30443 

0.95672 

0.07894 

-1.26770 

-1.19024 

6.00 

0.41543 

0.62092 

0.81671 

0.11663 

-1.25536 

-1 .18926 

5.00 

0.90382 

1.20144 

0.45889 

0.09916 

-1.23487 

-1.44020 

4.50 

1.24437 

0.96841 

1.19373 

-0.13653 

-1.20897 

-2.06755 

4.00 

1.64453 

-1.27575 

5.47300 

-0.84619 

-0.86801 

-  .30554 

3.75 

1.88349 

-1.13125 

6.45764 

-0.90367 

-0.78769 

.40946 

< 
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Table  VI.  Quantities  for  spin-orbit  matrices 


a  =  1/2 


2  + 

T 


V 


y/T  X 


/Fa 


-a 


fi  »  3/2 


Ga  A  =  .001255  au 
In  A  =  .00336 


.011835 
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Table  VII. 

Spin -orb it 

coefficients  for  the 

n  =  1/2  states 

R(a0) 

cr 

cn 

GaKr  3.75000 

.99955 

.02998 

4.00000 

.99960 

.02813 

4.50000 

.99961 

.02787 

5.00000 

.99927 

.03811 

6.00000 

.99641 

.08468 

7.00000 

.98199 

.18892 

8.00000 

.93379 

. 35782 

10.00000 

.82915 

.55902 

15.00000 

.81649 

.57736 

InKr 

3.75000 

.99704 

.07690 

4.00000 

.99738 

.07238 

4.50000 

.99742 

.07175 

5.00000 

.99534 

.09644 

6.00000 

.98063 

.19588 

7.00000 

.93712 

. 34900 

8.00000 

.87677 

.48091 

10.00000 

.82132 

.57047 

1 5 . 00000 

.81649 

.57736 

TIKr  3.75000 

.97505 

.22197 

4 . 00000 

.97731 

.21180 

4 . 50000 

.97763 

.21034 

5.00000 

.96502 

.26219 

6.00000 

.91737 

. 39803 

7.00000 

.86679 

.49867 

8.00000 

.83612 

. 54854 

10.00000 

.81788 

.57540 

15.00000 

.81649 

.57736 
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Table  VIII.  Rydberg-valence  transition  moments  in  GaKr 
(with  spin-orbit  corrections) 


III  1/2  -  I  1/2  HI  1/2  -  n  1/2 


III  1/2  -  I  3/2 


2  (X,Y) 


15.00  -0.75641  0.74574 
10.00  -0.72540  0.75618 
8.00  -0.44047  0.84363 
7.00  -0.22486  0.88026 
6.00  -0.10071  0.88449 
5.00  -0.05489  0.87255 
4.50  -0.05762  0.85454 
4.00  -0.00859  0.61353 


2  (X,Y)  (X,Y) 


-1.06969  -0.52733  -0.91335 
-1.07592  -0.50983  -0.91200 
-1.14948  -0.32327  -0.90345 
-1.16880  -0.16935  -0.89640 
-1.18499  -0.07517  -0.88767 
-1.43915  -0.03328  -0.87318 
-2.06674  -0.02383  -0.85487 
-0.30542  -0.01727  -0.61378 
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Table  IX.  Atomic  states  of  Ga  ,  In  ,  and  T1 

Excitation  Energy 

State  Ga  In  T1 


cm'^ 

eV 

cm”^ 

eV 

cm”  ^ 

eV 

2p 

1/2 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

3/2 

826.24 

.10241 

2212.56 

.274228 

7792.7 

.965840 

is 

i/2 

24788.58 

3.07234 

24372.87 

3.020814 

26477.5 

3.281665 

I.P. 

48380. 

5.9963 

46669.93 

5.784348 

49264.2 

6.105886 
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Table  X.  Rydberg-to-valence  transition  moments  In  model  TIKr 
(with  spin-orbit  corrections) 


R 

III  1/2 

Z 

-  I  1/2 
(X,Y) 

III  1/2 

Z 

-  II  1/2 
(X.Y) 

III  1/2 

(X.Y) 

15.00 

-0.75641 

0.74574 

-1.06969 

-0.52733 

-0.91335 

10.00 

-0.74665 

0.74590 

-1.06130 

-0.52476 

-0.91200 

8.00 

-0.67524 

0.75539 

-1.02925 

-0.49558 

-0.90345 

7.00 

-0.59354 

0.77699 

-1.03169 

-0.44701 

-0.89640 

6.00 

-0.47336 

0.81433 

-1.09099 

-0.35332 

-0.88767 

5.00 

-0.37761 

0.84264 

-1.38982 

-0.22894 

-0.87318 

4.50 

-0.43489 

0.83575 

-2.02130 

-0.17981 

-0.85487 

4.00 

-0.06471 

0.59985 

-0.29861 

-0.13000 

-0.61378 

I  3/? 
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Table  XI.  Horse  fitting  parameters : 

V(R)  -  V (**»)  =  Dg[u2  -  2u]  where  u  =  exp  [B(Rp  -  R)] 


Molecule 

State 

D 

e 

W 

«(•„ 

GaKr 

1 1/2 

.000772 

7.56 

.724 

1 3/2 

.00147 

7.19 

.72  2 

1 1 1/2 

.000625 

10.14 

.471 

1 1 1 1/2 

.00295 

6.26 

.853 

InKr 

1 1/2 

I3/28 

.000286 

7.86 

.817 

1 1 1  /2 
nIl/2a 

.000474 

10.24 

.483 

TlKr 

1 1/2 
h/2 

.00439 

8.52 

.633 

aSame  as  for 

1 1 1/2 

I 1 1 1/2* 

GaKr 

.000408 

10.12 

.502 
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Fig.  1.  Orbital  diagrams  for  the  electronic  states  of  Ga  +  Kr  and  Ga  +  Kr. 

The  lobes  and  circles  represent  the  in-plane  and  out-of-plane  d  orbitals; 
the  dashed  circle  denotes  the  Rydberg  orbital. 


THE  LOW-LYING  ELECTRONIC  STATES 
OF  GaKr  AND  GaKr  + 


RYDBERG 


Figure  2. 


THE  LOW-LYING  STATES  OF  GoKr  AND 


01 POLE  MOMENT  <0E8VES> 


(a) •'* 


38 


Figure  5. 


THE  LOW-LYING  STATES  OF  GaKr  AND 
GaKr+  WITH  SPIN-ORBIT  CORRECTIONS 
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Figure  6- 


DIPOLE  TRANSITION  MOMENTS  AMONG 
THE  LOW-LYING  STATES  OF  GaKr  WITH 
SPIN-ORBIT  CORRECTIONS 
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Figure  7. 


THE  LOW-LYING  STATES  OF  InKr  AND 
XnKr+  WITH  SPIN- ORBIT  CORRECTIONS 


R(A) 
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Figure  8. 


THE  LOW- LYING  STATES  OF  TIKr  AND 
TIKr+  WITH  SPIN- ORBIT  CORRECTIONS 
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Figure  9. 


DIPOLE  TRANSITION  MOMENTS  CONNECTING 
THE  LOW-LYING  STATES  OF  TIKr  WITH 
SPIN-ORBIT  CORRECTIONS 
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» 


Figure  11. 
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§3 


Application  of  the  Gordon-Kim  theory  to  the  group  I I IB-rare 
gas  systems. 

In  recent  years,  the  electron  gas  methods  developed  by  Gaydaenko 
and  Nikulin  [1],  and  Gordon  and  Kim  [2]  and  modified  by  Rae  [3]  and  Cohen 
and  Pack  [4]  have  proved  successful  in  calculating  the  interaction  energies 
of  pairs  of  closed  shell  atoms  or  ions  [1-5],  of  atom-molecule  systems[6-7] , 
of  molecule-molecule  systems  [8],  and,  recently,  of  systems  with  one  open 
shell  atom  [9].  These  interaction  energies  are  much  more  accurate  than  the 
simplicity  of  the  Gordon-Kim  (GK)  method  would  lead  one  to  expect. 

The  recent  success  of  the  GK  method  in  calculating  closed  shell  -  open 
shell  interaction  energies  has  prompted  us  to  use  this  method  to  calculate 
the  interaction  of  group  1 1  IB-noble  gas  pairs.  The  GK  theory  has  been 
most  successful  in  cases  where  the  interaction  is  non-covalent  and  where 
the  atoms  are  relatively  undistorted.  For  these  reasons,  systems  such  as 
GaKr  would  seem  to  be  ideal  candidates  for  a  GK  calculation. 

The  electron-gas  theory  is  briefly  reviewed  in  part  A,  and  our 
preliminary  results  on  GaKr  are  given  in  part  B. 

A.  The  electron  gas  method 

The  method  we  have  used  Is  that  of  Gordon  and  Kim  [2]  as  modified  by 
Cohen  and  Pack  [4].  A  more  detailed  description  of  the  theory  is  available 
in  these  two  papers.  Briefly  the  GK  theory  approximates  the  Intermolecular 


47 


potential  V(R)  at  the  distance  R  by 
V(R)  -  VGK(R)  =  V^p(R)  +  V^pr(R) 


where 


GK 

% 


=  vk  +  vc  + 


and  these  three  terms  represent  the  kinetic.  Coulomb  and  exchange  inter¬ 
action  energies,  respectively.  To  calculate  these  interaction  energies, 
the  electronic  charge  density  o  is  approximated  by  the  sum  of  the  atomic 
charge  densities, 


p  =  ca  +  Pb  • 

With  this  approximation,  the  Coulomb  interaction  can  be  calculated  directly, 
but  the  other  terms  are  all  estimated  by  the  formulas  for  the  energy  density 
of  a  uniform  electron  gas  [4].  An  additional  modification  [3,4]  is  made 
to  the  exchange  energy  to  avoid  self-exchange  contributions. 

8.  Calculations 

We  have  modified  the  molecule-molecule  GK  Interaction  program  of 
Parker,  Snow  and  Pack  [8],  which  allows  for  non-spherical  potentials,  to 
calculate  interaction  energies  for  closed  shell  -  open  shell  atomic  pairs. 
The  density  of  the  open  shell  atom.  In  this  case  Ga  ,  is  divided  into 
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the  spherically  symmetric  core  density  plus  the  valence  density.  For  Ga 

2  1 

in  the  ground  state  (4s  4p  ),  the  valence  density  is  constrained  to  be  in 
a  p  orbital  directed  along  the  internuclear  axis,  to  form  a  z  molecular 
state,  or  perpendicular  to  it  to  form  a  n  state. 

In  order  to  expedite  the  evaluation  of  integrals,  the  atomic  density 
is  expanded  in  a  set  of  Slater  basis  functions  with  the  coefficients  determined 
by  a  least-squares  fit  [4].  The  basis  set  expansions  for  the  density  are 
then  read  into  the  GK  program  which  calculates  the  interaction  energy 
by  three-dimensional  numerical  quadriture. 

We  have  obtained  numerically  tabulated,  relativistic  Hartree-Fock 
densities  for  the  group  1 1 1 B  and  rare  gas  atoms  from  Joseph  Mann  [10]. 

9  1  9  1 

The  densities  of  the  4s  4p  and  4s  5s  states  of  Ga  and  the  ground  state 
of  Kr  were  fit  with  small  sets  of  Slater  functions.  At  this  time  the 
fit  of  the  basis  set  expansion  is  not  very  good  (-v  10%). 


The  basis  set  expansions  obtained  from  this  fitting  procedure  were 
then  used  to  calculate  the  GK  interaction  energies.  Because  the  fitting 
procedure  does  not  normalize  the  density,  the  inaccuracy  of  the  present 
expansions  results  in  a  spurious  Coulomb  repulsion  between  the  two  atoms 
We  hope  that  improving  the  basis  set  fit  will  correct  this  error. 
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The  GK  potential  curves  are  compared  with  the  Cl  results  (see  Section  II) 

in  Figures  1-3.  The  energies  are  plotted  with  respect  to  the  asymptotic 

energy  of  each  state.  The  curve  labeled  GKR  includes  Rae's  correction  while 

the  GK  curve  is  the  unmodified  GK  theory.  The  GKR  curve  for  the  1^)+  state 

agrees  remarkably  well  with  the  Cl  result.  Unfortunately,  the  agreement  for 
2  2  + 

both  the  1  n  and  2  T.  states  are  not  as  good.  Specifically  the  depths  and 
positions  of  the  potential  wells  are  not  predicted  accurately.  We  are 
currently  working  on  improving  the  basis  set  fit  and  thereby  the  interaction 


energies . 
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§4 


4.1  Theoretical  basis 

A.  Introduction 

This  section  describes  the  theoretical  basis  of  our  calculations  using  the 
effective  potential  method.  For  an  excimer  system  AB,  where  A  is  a  closed-shell 
system,  most  of  the  states  of  interest  correspond  to  the  asymptotic  situation  where  R 
is  excited  but  A  is  in  its  ground  state.  Fundamentally,  what  the  electron  scattering 
theory  says  is  that  any  charged  particle  in  B  sees  a  potential  (e^),  due  to'  the 
closed-shell  system  A,  which  is  the  same  as  if  the  charged  particle  were  scattered 
off  of  A.  This  scattering  potential  is  corrected  for  the  fact  that  A  is  "de¬ 
polarised"  relative  to  the  scattering  problem,  due  to  the  presence  of  the  nucleus 
and  other  electrons  of  B.  The  unique  aspect  of  this  work,  which  distinguishes  it 
from  many  other  high-quality  attempts  to  produce  such  curves,  is  that  it  is  based 
on  the  model  interaction  potentials  and  response  functions  that  arise  out  of  many 
body  theory  (using  Schwinger  Functional  Derivatives).  The  use  of  this  type  of  many 
body  theory  has  characterized  the  work  of  this  group  over  the  last  few  years. 16 

The  effective  potential  method  will  clearly  give  incorrect  results  when  any  of 
•  three  situations  are  met:  First,  if  charge  transfer  states  mix  with  the  excited 
states  of  interest;  second,  when  states  go  asymptotically  to  two  excited  atoms; 
third,  when  covalent  bonding  is  important.  Only  for  charge  transfer  do  we  have  a 
possible  Class  I  solution.  For  the  moment  it  is  hoped  that  all  three  problem  situations 
will  be  avoided  for  the  lowest  excimer  excited  states  in  the  systems  of  interest  to 
this  proposal.  If  Zng  and  Mg^  are  any  guide  the  first  two  problems  occur  only  for 
the  higher  excited  states.17 

In  part  B  the  basic  equations  are  presented  and  possible  methods  for  calculating 
the  interaction  energy  are  discussed  in  part  C.  The  use  of  semi  empirical  forms  for 
the  effective  potentials  are  discussed  In  part  D. 
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8.  The  Basic  Equations 

Using  many-body  field  theoretic  methods  it  has  been  shown  that  the 
change  in  energy,  e  ,  resulting  from  the  addition  of  an  electron  to  a 
closed-shell  reference  system  (referred  to  here  as  A)  is  given  by  the  one- 
particle  Oyson  equation 

T(r)$(r)  +  f dr '  £A(r;r'  :e)4>(r)  «  e$(r)  (1) 

where  T  is  the  kinetic  energy  operator  and  <p  is  the  Dyson  amplitude 

with  r  and  r‘  being  space-spin  coordinates.  Thus,  the  problem  reduces 

to  an  effective  one-particle  problem  in  which  this  particle  experiences 

an  effective  potential,  eA  ,  which  represents  all  the  other  particles 

collectively,  taking  into  account  all  effects  such  as  polarization, 

correlation  and  exchange,  etc.  As  might  be  expected,  the  cost  of  this 

a 

simple  formulation  is  that  r  is  an  extremely  complicated  entity  which 
is  both  nonlocal  and  energy  dependent  and  which  cannot  rigorously  be 
brought  into  closed  form..  However,  it  has  been  possible  to  develop 
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excellent  closed  form  approximations  to  this  potential  which  are  based  on 
well  founded  physical  concepts.  Most  notable  among  these  is  the  Random 
Phase  Approximation  (RPA)  potential,  £Rp^  ,  which  has  been  very  successfully 
used  in  calculating  the  ionization  potentials,  excitation  energies, 
oscillator  strengths,  and  elastic  -scattering  phase  shifts  for  He  £133. 
iRpA  has  also  been  used  to  accurately  calculate  the  ionization  and  ex¬ 
citation  energies  of  Li  £193  Moreover,  it  has  been  shown  that  this  ab 
initio  potential  encompasses  other  phenomenologically  derived  semiempirica1 
potentials  which  have  been  used  by  other  workers  with  great  success  £203- 


It  has  also  been  recently  demonstrated  that  by  applying  the  same  many- 
body  techniques  to  the  problem  of  two  electrons  added  to  a  closed-shell 
reference  system  (A)  one  obtains  an  effective  two-particle  equation  for  - 
the  resulting  change  in  energy  £213- 


/dr'EA(ri;r':e))  +  +  /dr,  (dr2 'wV,  ,r2;r'r£)  ]*(r, 


e*(r,,r2)  , 


(2) 


where  in  addition  to  the  individual  one-particle  potentials  there  now 
appear*--  a  two-particle  effective  potential,  ,  which  represents  how  the 
presence  of  one  particle  affects  the  potential  seen  by  the  other  particle 
and  which  reflects  the  fact  that  system  A  can  act  as  a  dielectric  medium 

to  shield  the  Coulombic  interaction  between  two  charged  particles.  As 

a 

with  £  ,  this  potential  cannot  be  rigorously  brought  into  closed  form. 
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However,  an  excellent  ab  initio  closed  form  approximation  to  this  potential 
can  be  obtained  with  the  Random  Phase  Approximation,  WRpA  ,  which  is  completely 
compatible  with  the  similarly  obtained  one-particle  potential,  £RpA  • 


If  we  now  proceed  to  the  case  of  adding  m  electrons  to  our  closed- 
shell  reference  system  it  follows  by  induction  and  can  be  proven  formally 
that  the  change  in  energy  is  given  by 


dr 


*  /drj  Efl(r-  ;r!  :e)j  ♦  £.([7^  *  /drj  drj  wV/j  ;rj  ,rj  :t)) 
j  drj  irf,  UA(rj,rj,rk;rj,rj,r^:e)  +  ...  */drj  ...  drM/A(r, . .  -r^rj . .  .rj,:e)]  »(r, . . .  rn 


'  rr(rr..rm)  , 


(3) 


where  our  notation  for  three-particle  and  higher  potentials  is  obvious. 


In  the  above  we  only  considered  the  addition  of  electrons  to  our  closed- 
shell  reference  system  A.  However,  we  can  also  add  nuclei  as  well;  the  only 
difference  being  that  when  acting  on  a  nuclear  coordinate  all  potential 
terms  involving  electron  exchange  must  of  course  vanish.  Since  this  can 
be  trivially  accomplished  by  choosing  an  appropriate  (artificial)  nuclear 
spin  coordinate,  we  can  immediately  generalize  (3)  to  include  both  nuclei 
and  electrons  by  allowing  the  particles  to  carry  different  charges,  z  , 

(where  z  =  -1  for  an  electron)  to  give 
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m 


[£(T(r.) 
i-1  1 


/dr;  EA(rt;r;:e,2i))*f,  *  /drjdrj/fr,  ,rj;rj ,rj 


C,  1 


z-  z-)) 


•]  f (r. 


•rm} 


t,(rr  ••r») 


(4) 


We  of  course  cannot  solve  (4)  since  the  potentials  involved  cannot  be  written 
in  closed  form.  However,  we  can  replace  these  potentials  by  their  RPA 
approximates  which  are  in  closed  form.  Furthermore,  realizing  that  we  are 
deriving  a  theory  for  intermolecular  forces  which  is  essentially  perturbative 
in  nature,  we  will  now  assume  that  all  three-particle  and  higher  potentials 
can  be  neglected  so  that  our  equation  becomes  simply 


[£(T(r.)  - 
i  =  l  1 


/dr;  EA(r(ir;:e,Zi))*E  '  *  fadrAl'rjirrrj!*^  )KMr,..r 


'  e',(rr--r*)- 


(5) 


Since  our  interest  lies  in  the  calculation  of  potential  energy  surfaces, 
what  we  really  want  is  an  equation  for  the  change  in  energy  when  m& 

*  electrons  and  mn  nuclei  are  added  to  A  with  these  nuclei  held  fixed 

at  specific  points  in  space  (which  we  will  refer  to  collectively  as  R 
with  the  actual  spatial  coordinate  of  nuclei  i  being  R^  ).  That  is, 
we  want  to  be  able  to  separate  electronic  motion  from  nuclear  motion  so 
that  this  quantity, £(R)  ,  will  be  the  total  potential  experienced  by 
these  nuclei  and  will  satisfy  the  equation 

"n 

[£t(RJ  *£(R)]  X(R)  »  eX(R) 
i-1  1 


(6) 
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where  X(R)  is  a  function  of  the  nuclear  coordinates  only.  As  it 
stands,  (5)  precludes  such  a  separation  because  of  the  energy  dependence 
of  the  potentials  involved.  To  overcome  this  problem  we  will  now  assume 
that  for  those  solutions  to  (5)  we  seek  this  energy  dependence  is  not 
strong.  Furthermore,  we  will  also  assume  that  the  response  of  A  to  the 
added  particles  is  instantaneous.  In  this  way  we  can  replace  the  non- 
adiabatic  energy  dependent  potentials  in  (5)  with  their  hermitian  energy 
independent  adiabatic  approximates  such  as  those  given  in  ref. 21.  With 
these  substitutions  (5)  does  become  separable,  and  by  taking 


Hrv..rJ 


J 


(7) 


where  f(r^.. .r  )  is  a  function  of  only  the  electronic  coordinates 
e 

we  find  that  <£(R)  is  given  by 


<f(R)  -  EA<R)  -  £  ^(R^J+L  (TR-^lh-+  WA(R.,R  ;zvz  )),  (S) 

i=l  i>j  1  J  J 


A 

where  E  is  given  by  the  equation 


"e  mn 

CH ( ri -  -  ’rm^*R)  +  2^  fdr^A(r;r')  -  £  /drjWA(ri  ;r l;z j ) ) 


m„ 


*  /dridr !“*( ri * r j lr  1 r j  1  ]  »<rl— ■  EA(R)  »(r,...rBe)  (9a) 


where  H  is  the  usual  hamiltonian  for  me  electrons  In  the  field 

of  ro„  fixed  nuclei 
n 
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H(r,...r  ;R) 

e 


"le  mn  zj  me 

l/T(ri) '  £  W1  *  5  W 


m 


m. 


(9b) 


and  where  all  one-  and  two-particle  effective  potentials  are  taken  to  be 
hermitian  adiabatic  approximates  to  the  true  potentials.  Note  also  that 
all  potential  terms  involving  the  mn  nuclei  are  now  written  as  local 
quantities  thereby  taking  into  account  the  previously  mentioned  fact 
that  there  are  no  exchange  (i.e.  nonlocal)  terms  in  the  potentials  when 
nuclei  are  being  considered. 


If  we  now  collectively  refer  to  these  nr  electrons  and  m 

e  n 

nuclei  which  have  a  fixed  internuclear  geometry  as  being  system  B  , 
then  (f'(R)  'is  the  change  in  energy  resulting  from  the  creation  of  system 
B  in  the  vicinity  (as  measured  by  R  )  of  A  .  If  R  +  »  then 

which  is  just  the  energy  B  itself.  Therefore,  the  intermolecular 
potential  of  the  system  A-B  as  a  function  of  the  separation  between 
A  and  B  is 

vJ(R)  *  Ej(R)  -  E°  *  V*  nuc(R)  00.) 


where 


B.nuc 


00 


i’Rj4*j> 


i 


(10b) 
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and  where 


c°  »  /o  ,  v> B,B  1 

B  tv¥ 


00c) 


is  just  the  electronic  energy  of  isolated  system  B  and  where  Eg(R) 
is  given  by 


e  .  n 

[H(rv..r  ;R)  +  £  Jdr  1  (e^  ;rl )  -  £  wV .  ,R,  ;r:;Z?) ) 

e  i=l  J  j=l  J 

me 

*  £.  J<‘rjdr/(r).rjlr|rj)]  fj(rr..rm  )  =  E*(R)  »*<.•,. ..rm  )  (lOd) 


Therefore,  given  that  we  know  Eg  and  that  we  have  available  good 
closed  form  adiabatic  approximates  for  and  our  problem  reduces 
to  finding  the  solutions  (or  rather  a  particular  solution)  to  (lOd). 

In  the  following  section  we  will  look  at  some  ways  of  doing  this. 
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A 

C.  Determination  of  Eg 

a 

Our  equation  for  Eg  is  of  the  form 

(H  +  UA)  4^  =  EA  vA  (Ha) 

where 

a  me 

u  <ri"-'VR)  ■  £ 

me 

♦  E ;  /drJdrjwVj.rjirj.rj) 

* 

me  me 

*  £pV)  +  (llb) 

i=l  i>j  16  1  * 

where  we  note  that  fA  is  an  eigenfucntion  of  the  hermitian  'hainil tonian ' 

O 

(H  +  U^)  and  that  H  is  the  electronic  hamil tonian  for  isolated  system  B  . 

If  we  now  explicitly  assume  that  UA  is  small  compared  to  H  (this  of 

course  was  implicit  in  our  derivation  of  l)A  in  the  first  place)  then 

can  be  regarded  as  being  a  small  perturbation  on  H  .  As  a  consequence, 

the  solution,  ,  to  (11)  which  we  seek  should  resemble  the  electronic 

wave function  for  isolated  system  B  ,  T°  ,  and  this  wavef unction  should 

a 

therefore  provide  us  with  a  proper  starting  point  determining  h»b  and 

A 

o 

1.  First-Order  Perturbation  Treatment 


Projecting  (11a)  against  and  normalizing  tA  to  unity  gives 
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e£  »  <yJ|H  +  uV£>  ;  =  1  ,  (12a) 

whereas  if  fg  is  a  self-consistent  solution  for  isolated  system  B 
its  energy  is- 


-  <fg|H|fJ>  ;  <*X>  *  1 


Oi,„0 


(12b) 


If  we  assume  that  IP  is  a  quite  small  perturbation  to  H  then  we  can 

a 

solve  for  E"  using  standard  first-order  perturbation  theory  to  obtain 


'B  B 


E°  +  <T®|UA|f®> 


(13) 


A  ... 


so  that  our  first-order  perturbative  expression  for  Vg  is  simply 


VB 


yA 

B,nuc 


+  <f°lUA|y°> 


(14) 


This  of  course  is  equivalent  to  assuming  that  the  wavefunction  for  system 
B  remains  virtually  unchanged  when  in  the  presence  of  A  . 


In  terms  of  the  one-  and  two-particle  density  matrices,  p°  and  , 

o 
B 


associated  with  v?  ,  (14)  becomes 


^  *  VB,nuc  +  ,/ dr  P  (r  )p-|(rir‘)  +  /drldr2  q12^rl  ,r2^p2^rl  ,r2;rlrP- 

r  *r  r-r 

(15) 

If  Is  a  single  Slater  determinant  so  that 
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p°(rrr2;rj,r*)  =  1/2Cp®(r1  ;r*  )p°(r2;r' )  -  ;r*  )p°(r2;r' )]  (16a) 


o°(r;r') 


£  Mr)  Mr') 

1  =  1  1 


(16b) 


where  {*.}  are  the  n  one-electron  spin-orbitals  comprising  Y°  , 

1  rS 


then  Vg  takes  on  the  form 


^B.nuc  *  (,yhlPAl*l'  *  1/2  * 


where  the  two-electron  integrals  over  q*2  are  written  in  the  standard 
<  1 2 II  12  >  notation. 


If  y°  cannot  be  written  as  a  single  Slater  determinant  but  can 

D 

be  written  as  a  linear  combination  of  determinants  involving  n  one-electron 
spatial  orbitals  { <*>9 }  then  (12b)  can  always  be  written  as 

EB  =  -E.  Di  +  .  .Z.  .  Di.i  <*1*jJ7T2J*k*l>  {18) 

1  *J  ■  » J  »  K  *  1 

where  {D^,D^}  are  fixed  coefficients  which  depend  on  the  precise  form 
of  fg  and  on  orbital  overlaps.  Since  is  of  the  same  form  as  H 
and,  like  H  ,  is  spin  independent  (recall  that  isolated  system  A  must 


be  closed-shell)  we  can  immediately  write  our  first-order  approximation 
for  Vg  as 
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A  0 

=  V„  +  Z  O'? 
B.nuc  .  .  i 
* 


<»?lpA|4'°>  + 


i ,j ,k,l 


)k! 

ij 


.0.0, 


,A 

112! 


<<tT4>T  I  qV-  t 


vv 


(19) 


where  all  operators  and  integrals  are  now  explicitly  takent  to  involve 
only  spatial  coordinates.  Therefore,  we  can  determine  <vfg  1 UA ! 'f g>  by 

using  the  general  energy  expression  (18)  but  employing  one-  and  two-electron 
integrals  over  the  operators  pA  and  q^?  instead  of  the  usual  integrals 

_L 

over  h  and  r^  • 


We  could  continue  this  perturbation  treatment  by  going  on  to  determine 
second-order  and  higher  corrections.  However  we  will  stop  here  and  next 
consider  a  self-consistent  approach  instead. 


2.  Self-Consistent  Treatment 

A 

Let  us  now  assume  that  while  U  is  a  small  perturbation  to  H 

it  is  not  small  enough  to  justify  a  simple  first-order  perturbation 

A 

treatment.  That  is,  we  will  now  assume  that  fg  can  still  be  taken  to 
be  functionally  the  same  as  T®  but  because  of  the  presence  of  A  the 
spatial  orbitals  themselves  distort  away  from  {$?}  to  a  meaningful 
extent.  The  problem  is  then  to  determine  these  new  orbitals  {■5^}  . 

Since  fg  is  an  eigenfunction  of  electronic  motion  satisfying  ^11)  this 
can  be  done  varlatlonally.  That  is,  we  can  determine  these  new  distorted 

A 

orbitals  by  requiring  that  Eg  be  stationary  with  respect  to  changes  in 
these  orbitals.  If  for  simplicity  (but  not  necessity)  we  assume  that 
To  can  be  written  as  a  Hartree-Fock  type  wavefunction  involving  orthonormal 

D 
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spatial  orbitals  (or  orthonormal  spin-orbitals  for  the  case  of  an 
Unrestricted  Hartree-Fock  wavefunction)  then  (18)  takes  on  the  simple 
Roothaan  form 


Eu  =  2 

tB  *■ 


n  n 

Z  f •  <$?|h|4>?>  +  Z  [a. .  J  o  o  +  b.;K  o  o] 

i=l  1  11  i,j  ij  Vj  1J  Vj 


(20a) 


where 


J.o.o 
<t>  -  <t>  . 
i  J 


J 


<»0  0|_L 

1  Jr,. 


I  $V> 


(20b) 


K  o  o 

Vj 


.  K  <J> .  > 

i  1  <J> . 1  i 

J 


<♦.  K.  *  > 
J1  J 


A0A0|  I  I .0,0 

i  J  r-|2  J  1 


(20c) 


and  where  {f . ;a , . ;b^ i }  are  fixed  coefficients  (for  the  case  of  a  multi- 

configurational  wavefunction  these  coefficients  are  simply  related  to  the 

variationally  determined  configurational  coefficients)  which  depend  on 

the  precise  form  of  .  Since  fg  is  being  taken  to  have  the  same 

form  as  y°  and  because  of  the  similarities  between  H  and  UA  mentioned 

a 

above,  we  can  immediately  write  our  equation  for  Eg  as 


EA  =  2 

tg  ^ 


E  Vhlh#lv  ♦  I  [•u 

I  I  I  fj 


(21a) 


A  A 

where  the  operator  h  and  g^  are  given  by 


u  .  _A  ,  „A  1  .  A 

h  p  •*«  r£*"« 


(21b) 
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a 

Therefore,  }  and  therefore  Eg  can  be  determined  using  standard 

basis  set  expansion  (LCAO)  SCF  techniques.  The  only  difference  is  that 

instead  of  using  the  usual  one-  and  two-electron  integrals  we  must  use 

A  A 

integrals  over  the  operators  h  and  g^  instead.  However,  since 
these  integrals  serve  only  as  imput  this  difference  is  transparent  to 
whatever  available  SCF  procedure  we  employ. 

3.  Configuration  Interaction  Treatment 

For  the  sake  of  completeness  it  should  be  mentioned  that  since 
n 

can  be  variational ly  determined  there  is  no  need  to  stop  at  the  SCF 

level  and  we  could  solve  for  this  function  as  a  Configuration  Interaction 

(Cl)  problem.  The  only  difference  form  a  standard  Cl  calculation  is 

that  instead  of  using  the  usual  one-  and  two-electron  integrals  we  must 

A  A 

use  integrals  over  the  h  and  operators  defined  in  (21b).  It 

should  be  pointed  out  however  that  any  solution  for  Yg  which  differs 
significantly  from  Yg  implies  that  for  that  solution  UA  can  no  longer 
be  regarded  as  a  small  perturbation  to  H  and  in  such  a  case  the  validity 

A 

of  the  approximations  made  in  our  choice  of  U  would  become  subject 
to  question.  This  of  course  also  applies  in  our  SCF  treatment  as  well. 

D.  The  Potentials 

A  A 

The  one-  and  two-particle  effective  potentials  E  and  W 
appearing  in  our  final  equations  in  section  II  are  hermitian  adiabatic 
approximates  to  the  true  field  theoretic  potentials.  As  we  have  mentioned. 
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such  potentials  can  be  obtained  in  closed  ab  initio  form  using  many-body 
theory  within  the  framework  of  the  RPA  approximation  and  taking  the 
adiabatic  limits.  However,  while  these  potentials  are  tractable  they  are 
nonetheless  quite  complicated  and  their  use  would  entail  considerable 
computational  effort.  In  view  of  the  perturbative  nature  of  our  theory 
it  is  reasonable  to  expect  that  we  could  use  potentials  having  simpler 
forms.  Such  simpler  potentials  can  be  obtained  by  making  moment  expansions 
of  the  RPA  potentials  and  truncating  these  expansions  in  a  physically 
meaningful  manner  [21].  When  this  is  done,  the  resulting  potentials  can 
be  cast  in  forms  which  are  very  similar  to  phenomenologically  derived 
semi empirical  potentials  which  have  been  used  by  other  workers  with 
considerable  success  [223.  Therefore,  it  would  seem  that  the  use  of  com¬ 
plicated  aib  initio  potentials  is  not  warranted  (although  we  do  reserve- 
the  option  to  do  so)  and  that  we  can  take  our  potentials  to  have  semi- 
empirical  forms  similar  to  those  used  by  Dalgarno  and  by  Victor  [22], 
namely 


/dr'  A r;r') 


A  a 

“  17T  +  /dr’  EHF(r;r,)  "  ^7f4W6(klr|) 

aA 

-  — -9“ 6  WQ ( k | r- 1 )  +  (aQ  +  a1  |r|  +  a2 | r |2 )e~k I r I 
2 1  r  | 

(22a) 


and 


Jdrj  dr£  WA(r1  ,r2;rj  ,r£) 


A  (22b) 

*  ~T — 7T, — “2"  W3(k|r!  |  )W3(k|r2|  )Pn  (cos  r12) 

lrll  1 r2 I 

aA 

~lriT3Tr2l3  W4(klr1l)W4(k|r2[)P2(cos  y12) 
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where  we  have  chosen  our  coordinate  system  to  be  centered  on  A  which 

a 

for  simplicity  we  now  take  to  be  an  atom  having  a  nuclear  charge  of  z 

at  r 


and  where 


wn(x) 

y12 

P,(x) 

A 

ad 

A 

aq 

k 


n 

(1-e"x  )  is  a  cutoff  function 
angle  between  vectors  r-j  and  r ^ 
legendre  polynomial  of  the  i  degree 
dipole  polarizability  of  A 

approximate  quadrupole  polarizability  of  A  (adjustable) 
approximately  1/2  rQ  where  rQ  is  the  effective  radius 
of  A  (adjustable) 
adjustable  monopole  parameters 


and 


EA  (r;r' ) 
HF 


nA 

E 

i*l 


♦A  (r*)[2-Pr  r,]  ^(r1) 
- Ir^rT1 - 


is  the  static  Hartree-Fock  potential  of  A  with  being  the  nft 

spatial  Hartree-Fock  orbitals  for  the  electrons  in  A  .  P  ,  is  the 

r  ,r 

permutation  operator  if  r  is  an  electronic  coordinate  whereas  P  ,  3  0 

r,r 
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if  r  is  a  nuclear  coordinate.  Note  that  all  quantities  are  now  purely 
A  A 

spatial  and  that  z  and  W  are  therefore  explicitly  spin  independent. 

In  (22a)  the  first  term  is  simply  the  potential  due  to  the  nucleus 

of  A  and  the  second  term  is  the  static  Hartree-Fock  potential  for  the 

n 

electrons  in  A  occupying  the  spatial  orbitals  .  The  next  two 

terms  in  (22a)  are  asymptotically  correct  induced  dipole  and  quadrupole 
polarization  potentials  which  die  off  rapidly  at  short  distances  from  A  . 
The  final  term  is  an  induced  monopole  term  which  serves  as  a  short  range 
correction  potential.  The  terms  in  (22b)  describe  an  asymptotically 
correct  dielectric  potential  which  properly  cancels  out  one-particle 
induced  dipole  and  quadrupole  polarizations  of  A  due  to  two  particles 
of  the  same  charge  when  these  particles  are  on  opposite  sides  of  A  .  ' 
That  these  potentials  represent  a  significant  simplification  over  the 
ab  initio  potentials  is  clear  in  that  our  two-particle  potential  is 
strictly  local  and  the  only,  nonlocal  term  in  the  one-particle  potential 
is  simply  the  usual  Hartree-Fock  exchange  potential.  However,  despite 
their  simplicity,  potentials  such  as  these  have  been  used  very  success¬ 
fully  for  a  variety  of  problems  in  the  past  and  should  therefore  be  quite 
adequate  for  our  purposes. 

As  it  stands,  (22b)  is  concise  and  to  the  point.  However,  one 
important  property  is  obscured.  This  is  that  our  two-particle  potential 
can  be  written  in  terms  of  one-particle  operators  only.  To  see  this 
we  need  only  express  the  angle  in  terms  of  the  spherical  angles 
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for  each  vector.  Thus,  by  making  use  of  the  expansion 

Pz(cos  y12)  =  P,(cos  e,)P  (cos  e2)  +  2  z  fejj  P*  (cos  - 

m=  i v  7 

[cos  m«(i^  cos  m$2  +  sin  nty-j  sin  m$2] 

where  P^(x)  is  an  associated  Legendre  polynomial,  we  find 
can  be  rewritten  as 

g 

Jdr^  dr^  WA(r1  ,r2;rj  ,r^)  =  -  i  Q^r^Q^^) 

where 

Q](r)  ■  1  r*)~^W3(k) r )  )Pj  (cos  e) 

Qz(r)  =  >/aJ~|r|'2W3(kirj )  p|(cos  e)  cos  $ 

Q3(p)  -  |r|‘2W3(k|r| )  p](cos  0)  sin  4> 

Q4(p)  =  Irl'V^D  P2^cos 

Q 5(r)  =  |r|"3W4(k|r|)  pJ(cos  9)  cns  * 

Q6(r)  =  1/3 yja*  |r|'3Wr(k(r[)  P2(cos  9)  sin  * 

Q7(r)  =  i/12  yjai^  (r|-3W^(k|r  | )  P2(cos  0)  cos  2<t> 
%(r)  =  1  n2yja*  } r | ~3W1  ( k | r | )  P^(cos  9)  sin  24> 


^P^cos  a2) 
(24) 

that  (22b) 

(25a) 


(25b) 
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Consequently,  the  implementation  of  this  method  involves  the 
addition  of  only  one-electron  integrals  to  existing  SCF  or  Cl 

A 

programs.  In  addition,  if  I  is  replaced  by  a  pseudopotential, 
only  the  atomic  orbitals  of  B  will  be  used  to  calculate  the 


4. II  Details  of  the  Calculation 


In  developing  a  program  package  to  implement  the  effective  potential 
method,  our  goal  was  to  make  the  package  as  general  as  possible.  To  achieve 
this,  we  have  made  our  effective  potential  program  EFFPOT  compatible  with  the 
CVB  11^  SCF  and  Cl  programs.  This  program  package  is  a  gaussian  basis 
function,  molecular  SCF  and  Cl  routine  which  will  enable  the  present  version 
of  EFFPOT  to  be  extended  to  large  molecules.  The  version  of  GVB  II  that  we 
are  using  has  an  option  for  replacing  the  atomic  cores  by  a  pseudopotential; 
however,  the  form  of  this  built-in  pseudopotential  is  different  from  those 
being  developed  for  large  atoms^!  The  modifications  needed  to  make  use  of 
available  pseudopotentials  would  be  relatively  easy. 

The  GVB  II  progam  package  consists  of  a  number  of  sequential  program 
steps  which  are  outlined  below. 

1.  Integral  generation  a)  PLYLABS  -  creates  list  of  unique 

integral  labels 

b)  PLYINTS  -  evaluates  integrals 


2. 

Integral  preprocessing 

a)  PlYIJLK' 

convert  output  of  PLYINTS 

b)  PLYPAIR, 

to  form  needed  for  input 

to  GVB  II 

3. 

SCF 

GVB  II 

4. 

Cl 

a)  CORTRN 

transforms  to  molecular  basis 

b)  Cl  PROG 
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To  perform  an  effective  potential  calculation,  the  program  EFFPOT  is  inserted 
between  steps  2a  and  2b. 

FFFPOT  calculates  the  additonal  integrals  required,  <4>^  | |<j>  .>  and 
<4^  lQn  I  j>  >  anc*  modifies  the  integral  tape  so  that  the 

insertion  of  an  effictive  potential  is  transparent  to  GVB  II.  These  integrals 
have  not  been  previously  evaluated  for  gaussian  basis  sets,  and  because  of  the 

A 

cut-off  functions  in  i  and  Qn  and  the  angular  dependence  of  Qn,  the  evaluation 
of  these  integrals  is  rather  complicated,  especially  for  two-center  integrals 
containing  basis  functions  for  i  >  o. 

Although  it  is  relatively  easy  to  construct  an  algorithm  for  evaluating 
these  integrals,  the  algebraic  detail  rapidly  becomes  overwhelming.  To 

3 

alleviate  this  problem,  we  used  the  algebraic  programming  system  REDUCE  to 
produce  the  FORTRAN  code  required  to  evaluate  these  integrals.  An  explanation 
of  the  algorithm  used  to  derive  the  integral  expressions  and  an  example  of 
the  input  to  REDUCE  are  given  below  for  one  of  the  more  complicated  integrals. 
The  integral  we  shall  consider  is 

<  PxE  I  Q$  I  PxF>  03) 

v/here  the  superscripts  indicate  the  atomic  center  on  which  the  (basis  or 
potential)  function  is  centered,  Q4  is  defined  by  equation  10b  and  Px  indicates 
a  "Px"  gaussian  basis  function.  This  integral  can  be  written  explicitly  as 


<  P,r  I  Qu  I  p/>  = 


<*/3  fdrxee  ‘  r„3  L7i,  (kr0)F^  (cos  &a  )xF  e 


I. 


(14) 


Using  the  relation 


Pz  (  COS  ) 


-  75  - 


3  z5-r0a 


we  can  write 


I  -  i  Id r  xt e^r,'u,(l<r0)20l  r^xf  e'/?r/ 
"3f(dTXfe’<t  UZ,  O<r0 )  r03  x,  e  ^  ' 


2  I,*  I,. 


In  order  to  obtain  the  term  zg  ,  we  will  set  up  a  fake  gaussian  on  center  B, 
-cr|  5 

e  ,  and  move  the  W4(krB)/rg  term  to  a  fake  center,  D;  after  evaluating  the 
resulting  integral,  we  will  take  the  limits  D  -»■  B  and  c+  0.  Using  the 


relation 


Zj  e  cr°  =  eC<  +21  e~Cr° , 


we  obtain 

l»m  I  tw  _i.  — ! —  ■  - ■ 

I,  =  c-»o  D-*G  2  dE^F^ 


T7  ^3] 


where 


I,  -  ]  Jr  <f*r‘e  cr‘  r0"rLy  krj  eflr\ 


(i9)  ; 


We  can  also  write 


T  =  - T 

-LX  2M/3  df,3fT 
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where 

I^fdT.-^V^ri’U^kr,). 

We  can  then  apply  the  rule  for  combining  gaussians  on  two  centers. 


(22) 


where  Y  =  <A  + } 

REF  is  the  distance  between  centers  E  and  F,  and  center  A  is  defined 


(22a) 

by 


(22b) 


We  then  use  the  Fourier  Convolution  Theorem  to  obtain  the  following  expressions 
for  I3  and  I4: 

-  p  x  -  na 
-  y  “ff  X  ft 


J  =  „  r  nep  „  £  n*o 


where 


3  -  e  •  -  e  u  -  Is- 

(  -  Sr  * 

=  J  dz  e  r'/  (Ju  (k rft) 


(23) 


(24) 


(  dr  r* 


and 


-  sCfl  Q  *■ 

Xv  =  e 


(25) 


where 


*  (  Jr  e  ***  r0‘3  (JH  (kr6 ) 


(26) 


rk6[  Jr  rAL>,(kr){e“Y(r'P,,o)"-  e'y(r+R^j 
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It  these  equations  6  =  y  +  c  and  G  is  the  center  obtained  by  combining 

Gaussians  on  A  and  D.  The  derivatives  and  limits  in  equation  18  are  obtained 

analytically  by  REDUCE  and  the  result  is  that  equation  18  is  written  as  a 

sum  of  integrals  with  forms  similar  to  those  given  by  equations  24  and  26. 

These  basic  integrals  are  evaluated  numerically  and  summed  to  give  the  result 
E  B  p 

for  <PX  |Q4  |  Px  > .  To  obtain  expressions  for  this  integral  in  the  case 
where  any  of  centers  E,F,B  coincide,  REDUCE  is  used  to  obtain  the  appropriate 
limit  (E  -*■  F,  for  example).  The  REDUCE  commands  used  to  evaluate  all  Q4 
integrals  are  given  in  Appendix  B.  The  output  from  this  REDUCE  code  is 
rearranged,  using  additional  REDUCE  commands,  to  produce  FORTRAN  code  which 
can  be  inserted  directly  in  the  integral  evaluation  subroutines. 

The  present  versions  of  the  EFFPOT  integral  routines  are  restricted  to 
two  centers  and  s  or  p  basis  functions.  Our  tests  of  this  program  on  LiHe 
are  described  in  the  next  section. 

The  integration  itself  is  done  by  a  32  point  GauB  integration 
with  the  IBM  library  program  DQG  32.  The  other  FORTRAN-programs 
of  the  EFFPOT  routine,  developed  by  ourselves,  are  contained  in 
Appendix  C. 
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4. Ill  Test  Calculations  of  EFFPOT  on  LiHe 

In  this  section  we  present  the  results  of  the  first  stage  of  our  tests 
of  the  effective  potential  method.  We  have  chosen  LiHe  as  the  test  system, 
and  in  this  stage  of  the  tests  we  assume  that  the  parameters  of  our  semi- 
empirical  potential  are  known. 

Values  of  ad  and  can  easily  be  obtained;  however,  it  should  be 

pointed  out  that  in  a  semi-empirical  potential,  each  parameter  plays  a 

dual  role.  Besides  describing  the  physical  effect  to  which  it  most  obviously 

corresponds,  each  parameter  also  serves  to  correct  for  the  deficiencies  of 

the  semi-empirical  model.  Consequently  in  the  best  fit  for  the  potential 

parameters,  the  values  of  «d»  «q>  etc.,  should  not  be  expected  to  equal 

the  physical  quantities.  The  values  for  the  cut-off  functions  and  for  the 

short-range  part  of  the  potential  are  harder  to  obtain.  Our  initial  desire 

12 

was  to  use  the  parameters  obtained  by  Peach.  She  calculated  the  parameters 
for  a  model  potential  describing  one  electron  outside  of  He  by  fitting  the 
parameters  to  electron-He  atom  scattering  data.  Her  model  also  included  a 
pseudopotential  to  represent  the  Hartree-Fock  potential  of  He;  this  fact 
makes  it  impossible  to  separate  the  pseudopotential  from  the  short-range 
•  terms  needed  for  the  effective  potential  (see  equation  7).  Consequently, 

for  our  test  calculations,  we  use  Peach’s  values  for  cy,  Oq1 ,  8  and  her 
cut-off  functions.  In  addition,  we  have  guessed  a  value  for  a0  (equation  7), 
so  that  the  short-range  potential  has  only  one  term. 

The  exact  forms  of  the  potentials  used  in  our  test  calculations  are 
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| dr'  C*(  r;r': 

z)  = 

Z  z. 
r  + 

Z  /dr'  E *r(r;r‘) 

+  z3 

~Yr 

aoe 

jaUjt 
-  2  r* 

LJ2(f3r)  ~  *  rfc  UJ2  (p  'r ) 

(27) 

1  dr, ' dr \  r,: 

r  •  r'  f' 

.*2, ,2,)  = 

Z  r*  r/  (/?r;  )X//3Q  )  P  ( cos  4 ) 

-of'  2, 

Zz 

o 

Where  W  is  defined 

n 

by 

2  r,3 

r,3  X.^rJXjf/arJRfcos 

(28) 


and 


X„(») 


LJn  (x)  =  [  Xn(V)] 

r1  _  xn*' 

I  -  exp(-x)  L  ~0  (nTT)! 

w.o  w 


(29a) 


(29b) 


C<y  =  (  ~  (ojSjrt  )  . 

The  values  of  the  parameters  used  are  given  in  table  1  . 

In  tables  2  and  3  and  figures  1  and  2,  we  compare  the  effective  potential- 

1  3 

SCF  (EP-SCF)  calculations  with  SCF  and  restricted  Cl  calculations  for  the 

1 3 

same  basis  set.  The  basis  set  is  that  of  Krauss,  Maldonado  and  Wahl 

The  restricted  Cl  calculations  presented  here  have  the  same  asymptotic 

limit  as  the  EP-CI  calculations,  but  do  not  fully  allow  for  distortion  of  He 

at  intermediate  distances.  In  comparing  these  results,  it  should  be  noted 

that  the  EP-SCF  results  have  the  same  R  +  °°  limit  as  the  SCF  calculation, 

but  that  at  intermediate  distances,  the  EP-SCF  results  go  as  1/R^  while 

the  SCF  calculations  cannot  give  this  dependence.  Also,  the  full  Cl  calculations, 

which  allow  for  more  He  polarization,  should  give  more  attractive  interaction 


Parameters  for  the  potential  defined  by  equations  27  to  29,  for  He' 


Krauss,  Maldonado,  and  Wahl 
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Krauss,  Maldonado,  and  Wahl 


Table  3.  LiHe  interaction  enrgies  (in  Hartrees 
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Figure  Captions 

1.  Interaction  energy  for  LiHe  2l  state. 

©  connected  by  solid  line;  EP-SCF  results. 

1 3 

H  SCF  results  of  Krauss,  Maldonado,  and  Wahl 

•  restricted  Cl  calculation. 

2.  Interaction  energy  for  LiHe  2n  state. 

Q  connected  by  solid  line;  EP-SCF  results. 

n  1 3 

LI  SCF  results  of  Krauss,  Maldonado,  and  Wahl 

•  restricted  Cl  calculation. 


Interaction  Energy  *103  (Harcrees) 
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4 . IVa  Test  calculations  of  the  fit  program 

In  fig.  2  in  the  last  paragraph  it  has  become  obvious  that  the 
initial  values  for  the  parameters  in  the  effective  potential  (the 
ones  given  in  table  1  in  section  III,  let  us  call  this  set  of  para¬ 
meters  in  the  following)  do  not  give  a  sufficiently  good  agreement 
between  the  calculated  molecular  potential  curve  and  the  curve  of 
our  Cl  calculation  or  the  results  of  Krauss  et  al.  One  method  to 
find  a  better  fit  would  be  to  repeat  the  calculations  for  a  large 
number  of  parameter  values  and  so  to  calculate  a  potential  curve 
for  each  point  on  a  dense  grid  in  the  10-dimensional  parameter  space. 
This  would  be  the  safest  method  to  come  close  to  the  best  possible 
fit  but  it  would  consume  an  enormous  amount  of  computer  time  and 
therefore  it  would  be  contrary  to  our  initial  goal.  A  faster  alter¬ 
native  is,  to  apply  a  fit  program  which  uses  as  input  the  results 
of  section  III  and  tries  to  improve  the  molecular  potential  curve  by 
some  kind  of  first  order  perturbation  calculation. 

In  order  to  explain  what  kind  of  fitting  method  has  been  used 
here  let  us  for  this  section  introduce  the  following  notation,  p 
is  a  set  of  lO  parameters  which  go  into  the  effective  potential, 

W(p,R)  is  the  complete  effective  potential  for  the  set  p  of  parameters 
and  for  the  nuclear  distance  R,  Y (p, R)  is  the  exact  electronic  SCF 
wavefunction  for  the  nuclear  distance  R  and  for  the  effective 
potential  W(p,R)  calculated  according  to  the  method  given  in 
section  III.  V(p,R)  is  the  molecular  potential  calculated  according 
to  the  method  given  in  section  III  with  the  effective  potential 
W (p , R)  and  at  a  nuclear  distance  R. 

The  idea  of  our  first  order  perturbation  approach  is  to  calculate 


first  the  difference  between  the  effective  potentials  for  two  sets 
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of  parameters  pQ  and  p  ,  then  calculate  the  expectation  value  of 
this  difference  potential  with  the  exact  wavefunction  belonging  to 
parameter  set  p,  and  interpret  this  expectation  value  as  difference 
in  the  molecular  potential.  In  formula 

V(p  ,  R)  =  V(Pq,R)  +  <iMpo,R)  [W(p,R)  -  W(Pq,R)  |4'(po>R)  >  (30) 

as  justification  for  the  application  of  a  first  order  perturbation 

method  we  have  first  to  show  that  the  results  of  the  exact  calculation 

are  highly  linear  at  variation  of  a  parameter  in  the  effective 

potential.  In  fig.  3  the  electronic  energy  has  been  calculated  with 

the  method  given  in  section  III  for  various  values  of  the  dipole 

parameter  and  always  the  same  values  for  the  other  parameters  at 

the  same  nuclear  distance  of  R  =  3.5.  It  can  be  seen  that  the 

dependence  of  the  energy  on  the  parameter  is  highly  linear.  Similar 

results  hold  for  the  other  parameters.  These  calculations  have  also 

shown  that  the  wavefunction  does  not  change  significantly  under  a 

change  of  a  parameter.  Therefore  it  seems  that  a  first  order  perturbation 

approach  to  the  fitting  problem  is  justified.  One  of  the  inputs  in 

any  fitting  program  must  be  a  value  to  which  we  want  to  come  close. 

In  the  case  here,  the  assumed  true  molecular  potential  curve  for  the 
2 

He-Li  tt  state  is  the  one  given  by  Krauss  et  al.  We  denote  it  in 
the  following  by  VK(R).  In  addition,  it  is  necessary  to  define  the 
quantity,  a  measure  of  the  error,  which  should  be  minimized. 

We  have  chosen  as  measure  of  the  error  of  the  molecular  potential 
curve  V(p,R)  the  following  quantity 


Dtp)  =1  [ V  ( R .  ) 

1  K  l 


V(p,R. ) ]2/V^(Ri) . 


(31  ) 
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The  sum  i  goes  over  all  those  nuclear  distances  at  which  we  try  to 
improve  the  function  V(p,R)  at  the  moment. 

This  fitting  for  the  relative  error  is  a  disadvantage  in  so  far, 
that  those  nuclear  distances  at  which  VR(R^)  is  small  have  a 
higher  weight  than  those  distances  at  which  V  (R  )  is  big.  Therefore 
the  computer  tries  to  fit  V(p,Ri)  with  the  higher  priority  near 
zeros  of  • 

The  fitting  program  and  all  subroutines  which  belong  to  it  and 
all  those  programs  which  are  necessary  in  addition  to  the  previous 
programs  are  shown  in  Appendix  D. 

Our  experience  running  this  fitting  program  can  be  summarized 
as  follows:  The  fitting  program  is  started  with  a  particular  set 
p.^  of  10  parameters  and  in  its  Oth  step  it  calculates  the  quantity 
V{pi#R)  according  to  (30).  For  this  initial  step  it  takes  about 
6  minutes  times  the  number  of  points,  on  which  the  curve  is  fitted, 
to  calculate  the  integrals.  Such  a  long  time  is  necessary  because 
up  to  now  there  are  no  optimised  standard  programs  of  quantum 
chemistry  available  for  these  types  of  integrals.  It  takes  a  long 
time  to  transform  the  expressions  of  the  functions  to  be  integrated 
into  such  a  form  that  a  library  integration  program  can  be  applied. 
(The  32  point  Gaussian  integration  program  DQG32  by  IBM  ultimately 
has  been  used.  If  in  a  next  step  one  of  the  parameters  a^,  a^,  6, 
aQ,  a^ ,  a^,  a^  is  changed  then  it  is  not  necessary  to  evaluate  the 
integrals  again  because  these  7  parameters  do  not  change  the  shape 
of  W(p,R).  They  go  linearly  into  W(p,R)  and  also  the  corresponding 
integrals  are  linear  in  these  parameters.  Therefore  it  takes  only 
a  few  seconds  to  calculate  a  new  V(p^,R)  when  p^  differs  from  p.^ 
only  in  a  parameter  which  goes  into  W(p,R)  linearly.  The  situation 


is  completely  different  if  one  of  the  cut-off  parameters  is  changed. 
Then  the  shape  of  the  function  W(p,R)  is  different  and  it  is  necessary 
to  calculate  from  the  start  the  corresponding  integrals  again.  This 
takes  at  least  several  minutes  time.  Hence  it  is  quite  time  consuming 
to  fit  cut-off  parameters  by  this  fitting  program.  Therefore  we  have 
done  the  fit  of  cut-off  parameters  partly  by  hand. 

In  Fig.  4  the  results  of  the  fit  can  be  seen  for  two  sets  of 

parameters  which  gave  the  best  fits  we  could  obtain  for  the  molecular 

2 

potential  curve  of  the  tt  state  of  HeLi.  For  comparison  we  show  in 

addition  the  curve  calculated  by  Krauss  et  al.  as  solid  line.  This 

curve  agrees  quite  well  with  our  Cl  calculations  (see  Fig.  2  in 

section  III) .  The  result  of  the  SCF  calculations  from  section  III 

is  shown  as  dotted  line  also  for  comparison.  This  line  is  the  molecular 

potential  curve  which  belongs  to  the  parameter  set  pQ  and  it  is  the 

Oth  order  result  of  the  fitting  procedure  in  this  section.  The  broken 

line  and  the  chain  line  show  the  1st  order  results  for  two  different 

sets  of  parameters.  It  is  obvious  that  the  deviation  between  these 

first  order  curves  and  the  Krauss  et  al.  curve  is  a  lot  less  than  the 

deviation  between  the  Oth  order  curve  and  the  Krauss  et  al.  curve. 

The  worst  error  of  our  1st  order  curves  is  in  the  region  near 

R=5  to  R=6 .  In  this  region  it  has  not  been  possible  to  get  a  curve 

which  is  steep  enough.  All  considered  1st  order  curves  with  a 

reasonable  over  all  shape  have  been  too  horizontal  in  this  interval. 

Up  to  now  we  did  not  yet  have  time  to  make  calculations  for  molecular 

2 

potential  curves  of  other  states  of  LiHe  than  the  tt  state.  Before 
we  do  this,  it  would  be  useful  to  try  out  a  new  fitting  strategy 
described  in  the  next  section. 
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T7b.  A  strategy  to  achieve  a  better  fit  to  the  tt  state  of  LiHe 

As  was  seen  in  IVa  we  did  not  get  an  ideal  fit  in  the  R=5  region 
as  we  did  not  seem  to  get  the  fitted  curve  to  be  steep  anough.  It 
is  to  be  noted  that  so  far  we  have  fitted  to  Krauss  et  al .  LiHe, 
it  state.  It  is  also  to  be  noted  that  we  have  not  exhausted,  because 
of  lack  of  computer  time,  the  fit  routine  and  tried  to  guess  the 
best  parameter  regions.  In  this  subsection  we  suggest  a  new  strategy 
which  we  plan  to  carry  out  in  the  fall  of  1980  when  U.S.C.  will 
replenish  our  computer  budget.  The  strategy  is  based  on  the 
comparison  of  the  shape  of  the  ground  LiHe+  potential  curve  (which 
goes  asymptotically  to  Li+  +He)  to  that  of  the  SCFtrLiHe  curve  (which 
dissociates  properly)  as  given  in  the  figures  in  Krauss  et  al.  In 
this  we  see  the  striking  fact  that  the  two  curves  are  parallel  in 
the  R=5  to  6  region  and  almost  parallel  at  short  range;  the  LiHe  has 
a  slightly  smaller  Re.  This  suggests  the  much  of  the  curve  in  tt  LiHe 

is  like  the  X1LLiHe+  and  this  is  clearly  charge-induced  dipole  at 

-4  + 

large  distance  (R>4.5)  as  it  goes  at  R  and  Li  -  Helium 

repulsion  at  short  distance.  Now  the  former  effect  is  in  our  model 

4 

due  to  the  a/r  term  and  the  latter  in  the  static  and  monopole  term, 
sir 

(in  which  e  ((a  chosen  so  that  the  term  effects  the  very  short 
range  part  of  the  curve  only))  is  ‘always  automatically  fit  with  a 
negative  coefficient) .  This  suggests  that  as  a  first  run  at  fitting 
the  tt  ,  that  we  work  with  the  ionic  state  (a  simpler  problem)  and 
set  our  quadrupole  term  to  zero  and  perhaps  again  start  with  the 
physical  value  of  the  polarizability  for  Helium.  This  simpler  system 
should  be  able  to  be  fit.  The  ionic  system  then  determines  the  a, 
its  short  range  cut-off  parameters  and  will  give  a  first  guess  of 
the  monopole  term. 
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Now  we  will  return  to  the  tt  L^-He  fit,  and  retain  the  LiHe+ 

parameters  as  a  start.  We  will  then  work  at  large  R,  R>7  only,  fit 

the  quadrupole  term  and  its  cut-off  parameters.  This  is  suggested 
2  1 

again  by  the  r.  state  to  X  I  ionic  state  comparison  figure  in  Krauss 
where  we  see  that  at  large  R  the  LiHe  system  deviates  from  the  ionic 
system  and  is  less  attractive,  i.e.  the  effect  of  a  R^  comes  in. 

This  effect  can  only  be  found  in  the  it  state  as  at  large  R,  the 
Li  P  state  can  be  reviewed  roughly  as 


i.e.  Helium  sees  a  {  +  }  charge  distribution  and  the  L  tt  orbit  feels 
the  Helium  quadrupole  moment.  The  next  step  will  be  to  now  fit  again 
at  small  R,  the  monopole  term,  starting  with  our  ionic  result  •  The 
small  shift  must  be  due  to  a  chemical  effect  which  can  be  described  as 

— 2fTT  Lr 
/-  2s 

H e  _ ' 

We  should  be  able  to  account  for  this  with  our  monopole  parameter. 

Hopefully  this  new  strategy  will  yield  a  better  fit  with  little 
computing.  The  effective  potential  will  then  be  used  to  construct 
the  z  state. 


/ 

Up 
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Failure  in  fitting  by  the  procedure  of  this  subsection  or 
in  the  construction  of  the  I  state  will  forebode  badly  for  the 
ultimate  success  of  this  theory. 
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Figure  Captions 


Figure  3 


2 

Plot  of  the  electronic  energy  for  the  He-Li  tt  state  calculated 
according  to  the  method  given  in  section  III  for  various  values  of 
the  paramter  and  otherwise  always  the  parameters  given  in  table 
1  in  section  III.  The  3Uclear  distance  is  always  R=3.5. 


Figure  4 : 

2 

Plot  of  potential  curves  for  the  He-Li  tt  state  calculated  with 
various  parameters  for  the  effective  potential  and  the  curve  given 
by  Krauss  et  al .  for  comparison. 

Solid  line:  result  of  Krauss  et  al . 

Dotted  line;  SCF  calculation  from  section  III  with  parameter  set  pc 
Broken  line:  calculation  with  the  method  given  in  section  IV  and 
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We  report  ab  initio  configuration  interaction  calcuUtioni  on  the  SUMS  of  the  gallium  krypton  (GaKr) 
molecule  arising  from  the  Ga(I/>,,j  j/i,1S1,i)  +  Kr(‘S0)  and  Ga*(‘S0)  +  Kr(‘S0)  separated  atom  limits.  The 
potential  energy  curves  for  the  states  arising  from  the  Ga(*P|;,j;i)  limits,  the  /  1/2,  II  t/2.  and  /  3/2 
states,  are  found  to  be  repulsive.  The  potential  energy  curves  for  the  states  arising  from  the  Ga(JS|,})  and 
Ga*('S0)  limits,  the  III  1/2  and  1  0  states,  are  both  found  to  be  weakly  bound;  D,  (III  1/2)  =  0.047 
eV  and  0(1  0)  =  0.24  eV.  The  potential  energy  curves  and  transition  moments  obtained  in  the  GaKr 
calculations  have  been  used  to  simulate  the  curves  and  moments  for  InKr  and  TTKr  Using  this  data  the 
absorption  and  emission  coefficients  of  all  three  molecules  have  been  calculated. 


I.  INTRODUCTION 

The  Group  IHA -rare  gas  exclmer  molecules  are  con¬ 
sidered  to  be  attractive  candidates  for  developing  an  ef¬ 
ficient,  high  power,  tunable  laser  in  the  visible  region 
of  the  spectrum.1  The  transition  under  consideration  is 
a  perturbed  (it+  U*SI/2-«,J>1/1,J/l  transition  of  the  Group 
niA  atom.  While  current  experimental  studies  have 
concentrated  on  the  thallium -rare  gas  'Sterns,  especial¬ 
ly  TIXe ,  as  the  most  promising  candidates ,  Gallagher* 
has  recently  raised  the  possibility  of  using  the  gallium- 
rare  gas  systems  by  obtaining  gallium  atoms  from  the 
dissociation  of  Gal, . 

To  provide  further  information  on  the  nature  of  the 
exclmer  states  involved  in  these  studies ,  we  report  here 
ab  initio  configuration  interaction  calculations  on  a  pro¬ 
totype  Group  IIIA -rare  gas  diatomic  molecule,  GaKr. 

We  present  the  potential  energy  curves  for  ail  of  the 
states  arising  from  the  neutral  Ga(*Pl/ltl/J,  *S1/a) 

+  Kr(‘5^)  and  ionic  GaVS^  +  Kr^S,)  separated  atom  lim¬ 
its  and  the  dipole  transition  moments  radlatively  cou¬ 
pling  the  states .  We  then  use  the  computed  potential 
curves  and  transition  moments  for  GaKr,  along  with  the 
experimental  atomic  spin-orbit  coupling  constants,  to 
model  the  curves  and  moments  for  InKr  and  TlKr.  With 
this  data  we  calculate  the  emission  and  absorption  coef¬ 
ficients  for  all  three  systems  using  the  classical  tech¬ 
nique  developed  by  Gallagher  and  co-workers.1* 

II.  ELEMENTARY  THEORETICAL  CONSIDERATIONS 

Let  us  first  consider  the  description  of  the  states  of 
the  Group  IllA-rare  gas  molecules  without  spin-orbit 
corrections.  The  valence  Ga(*P) ♦  Kr(‘S)  separated  atom 
limit  gives  rise  to  a  *n  state  and  a  *2*  state,  the  Ryd¬ 
berg  Ga(,S)  +  Kr(‘S)  limit  gives  rise  to  another  *2‘  state 
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and  the  Ionic  Ga*(,S)  +  Kr(‘S)  Limit  gives  rise  to  a  '2* 
state .  We  will  label  the  states  which  arise  from  the 
valence  limit  the  1  *n  and  1  *2*  states ,  that  from  the 
Rydberg  limit  the  2*2*  state  and  that  from  the  ionic  lim¬ 
it  the  1  ‘2*  state. 

The  orbital  diagrams  for  the  1  2D ,  X  *2*,  and  1  '2* states 
of  a  Group  IIIA -rare  gas  molecule  are  given  in  Fig.  1. 
From  these  diagrams  it  is  evident  that  none  of  the  re¬ 
sulting  potential  energy  curves  are  expected  to  be  chemi¬ 
cally  bound.  In  the  ionic  state,  however,  weak  binding 
can  result  from  the  charge -induced  dipole  interaction  of 
the  Group  EDA  ions  and  the  rare  gas  atoms.  In  fact,  this 
Interaction  might  also  be  expected  to  give  rise  to  a  weak 
binding  in  the  2*2*  state.  Since  the  Rydberg  orbital  of 
the  Group  IIIA  atom  is  diffuse,  the  approaching  rare  gas 
atom  can  easily  polarise  the  Rydberg  orbital  out  of  the 


THE  LOW-LYING  ELECTRONIC  STATES 
OF  GaKr  AND  GoKr* 


Go  Kr 

4=4° 


FIG.  1.  Orbital  diagrams  for  the  low-lying  electronic  states 
of  GaKr  sad  GaKr*.  The  two  lotted  figures  represent  4 p  orbi¬ 
tals  in  lbs  plane  of  the  paper;  the  circle  represents  a  ip  orbi¬ 
tal  perpendicular  to  the  plane  of  the  paper;  the  Ss  Rydberg 
orbital  Is  rsprsssaisil  by  a  large  dashed  circle. 
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interatomic  region,  thus  partially  unshielding  the  Group 
UIA  ionic  core.  The  binding  in  the  2*E*  state  should, 
of  course,  be  substantially  less  than  that  in  the  1  'e' 
state. 

For  the  valence  states  the  repulsiveness  of  the  curves 
should  be  roughly  proportional  to  the  number  of  po  elec¬ 
trons  since  the  overlap  of  the  two  atomic  charge  distri¬ 
butions  is  dominated  by  the  overlap  of  the  po  orbitals. 
Thus,  the  potential  curve  (or  the  1  *11  state  with  two  po 
electrons  should  be  less  repulsive  than  that  for  the  1 
state  with  three  po  electrons. 

As  we  shall  later  see,  at  short  internuclear  separa¬ 
tions  the  energies  of  the  1*E*  and  2*E'  states  are  nearly 
equal  and  there  is  a  strong  interaction  between  the  two 
states.  This  mixing  is  strong  in  spite  of  the  fact  that 
one  state  is  a  valence  state  and  the  other  a  Rydberg 
state  because  the  two  zero-order  configurations  (see  the 
next  section)  differ  by  only  a  single  excitation. 

A  complete  treatment  of  the  electronic  states  of  the 
Group  IUA-rare  gas  molecules  must  include  the  effects 
of  spin-orbit  coupling.  Using  the  simple  model  devel¬ 
oped  in  earlier  papers  for  including  the  effects  of  spin- 
orbit  coupling  in  diatomic  molecules,1  only  the  1*11  and 
l*E*  states  are  affected  by  the  spin-orbit  interaction, 
and  the  coupling  arises  solely  from  the  spin-orbit  inter¬ 
action  in  theGroup  OlA  atom.  For  gallium  the  *P,/2- 
tP]/l  splitting  is  only  0.10  eV*;  however,  by  thallium 
this  splitting  has  increased  to  0.97  eV.’ 

In  a  molecule  only  the  z  component  of  the  total  angular 
momentum , 

O-A  +  S,  , 

is  a  good  quantum  number.  Thus,  from  the  Ga(*P 1/JtJ/,) 
+  Kr(,S#)  separated  atom  limits  we  obtain  two  0=  1/2 
states  and  one  ft  =  3/2  state  (the  /  1/2,  11  1/2,  and  f  3/2 
states),  from  the  Ga(*S,„)  +  Kr('S,)  limit  another  0=1/2 
state  (the  111  1/2  state),  and  from  the  Ga*(tS,)  +  Kr(\S,) 
limit  an  0  =  0  state  (the  I  0  state).  At  large  internuclear 
distances  the  f  3/2  and  II  1/2  potential  curves  are  degen¬ 
erate  and  are  separated  from  the  f  1/2  (ground  state) 
curves  by  the  group  IIIA’P, n-*P\n  splitting.  As  the 
distance  decreases  the  1  *fl  and  1  *£*  curves  rapidly 
separate,  the  separation  Increasing  approximately  ex¬ 
ponentially.  As  a  result  at  short  distance  the  atomic 
aptn-orbit  coupling  is  nearly  quenched  and  the  1  1/2  and 
11  3/2  curves  approximate  the  1  *11  curve  and  the  U  1/2 
curve  approximates  the  1  *E*  curve.  The  potential 
curves  of  the  2*E*  (111  1/2)  and  1  'e*  (f  0)  states  are  un¬ 
affected  by  the  Inclusion  of  spin-orbit  coupling. 

III.  OETAILS  OF  THE  CALCULATIONS 

A.  Basis  sets 

The  calculations  employed  an  atomic  (MsllpSrf)  primi¬ 
tive  Gaussian  basis  eet  lor  both  the  galliam  and  krypton 
atoms.’  The  core  orbitala  (1*  ,2*,3*,2p,3p,  and  3d) 
were  contracted  to  a  single  function  while  the  valence 
orbitala  (4s  and  4F)  wera  contracted  to  two  functions 
using  the  general  contraction  method  of  Raffenetti*  (aee 
also  Rsi.  7).  These  atomic  basis  sets  were  segmented 


with  a  single  set  of  3d  functions,  f„.  -  0. 16  and  tKr 
■  0.35,  to  describe  molecular  polarization  effects.  The 
values  of  these  exponents  were  based  on  calculations  on 
other  excimer  laser  systems.  *•* 

Since  the  transition  of  interest  involves  the  first  Ryd¬ 
berg  state  of  the  gallium  atom,  the  gallium  basis  set  was 
augmented  with  two  diffuse  $  functions,  t,  =  0.026  and 
0.011,  and  a  diffuse  p  function,  £,  =  0.010,  to  describe 
this  state.  The  exponents  of  these  functions  were  ob¬ 
tained  from  Hartree-Fock  (HF)  calculations  on  the  *S 
and  *P  Rydberg  states  of  the  gallium  atom. 

The  final  basis  sets  thus  consisted  of  a  (16sl2p6rf) 
primitive  set  contracted  to[7s5/>2d]  for  gallium  and  a 
(14sil/>6d)  primitive  set  contracted  to  [5*4/>2dJ  for 
krypton. 

B.  Calculations  on  the  gallium  atom 

The  reference  orbitals  for  the  configuration  interac¬ 
tion  (Cl)  calculations  on  the  gallium  atom  were  obtained 
from  a  Hartree-Fock  (HF)  calculation  of  the  *P  state  of 
the  atom  with  the  configuration 

*P:  . .  •  4s*4p  .  (1) 

(The  4 P  set  includes  all  three  cartesian  components  and 

all  are  equivalent.)  The  Rydberg  5s  orbital  was  ob¬ 
tained  from  an  IVO  calculation’  using  the  4s*  core.  In 
terms  of  these  orbitals  the  reference  configurations  for 
the  other  states  of  interest  are 

*S:  •  •  •  4s*5*  (2) 

‘S:  4s*  .  (3) 

Both  polarization  Cl  (POL-CI)10  and  full  Cl  calcula¬ 
tions  have  been  carried  out  on  the  above  states  of  the 
gallium  atom.  The  POL-CI  calculations  include  all 
■Ingle  and  double  excitations  relative  to  configurations 
(l)-(3)  with  the  restrictions  that 

(1)  all  core  orbitals  (ls-3s,  2p-3p,  and  3d)  remain 
doubly  occupied  and 

(2a)  no  more  than  one  electron  occupy  the  Rydberg  5s 
and  virtual  orbitals  (*P  and ’s  states)  or 

(2b)  no  more  than  one  Electron  occupy  the  Rydberg  5s 
orbital  and  no  more  than  one  electron  occupy  the  virtual 
orbitals  (*S  state). 

The  POL-CI  calculations  involve  22  space  and  33  space/ 
spin  configurations  for  the  *P  state,  27  space  and  43 
space/spin  configurations  for  the  *5  state,  and  15  space 
and  15  space/spin  configurations  for  the  'S  state. 11  Note 
that  the  POL -Cl  wavefunctlon  accounts  for  the  4s*- 4p* 
•ear  -degeneracy  effect. 

The  less  restrictive  condition,  (2b),  for  the  *$  state  is 
necessary  to  obtain  a  comparable  description  of  this 
state.  Only  by  relaxing  condition  (2a)  can  configurations 

such  as 

•  *  •  l 

be  included  in  the  POL-CI  wavefunctlon  of  the  *S  state. 
These  configurations  account  for  the  first -order  differ¬ 
ence  between  the  41*  orbital  obtained  from  the  calcula- 
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ticns  on  the  *P  state  and  the  p  orbital  needed  to  describe 
the  angular  correlation  of  the  4s*  pair  (4s-4p  near  de¬ 
generacy)  and  are  quite  important.  Comparable  config¬ 
urations  are  automatically  included  in  the  POL-C1  wave- 
functions  of  the  *P  and  'S  states . 

The  full  Cl  calculations  include  all  excitations  (single - 
triple)  relative  to  the  configurations  (1)— (3)  with  only 
restriction  (1)  above.  The  full  Cl  calculations  include 
157  space  and  263  space/spin  configurations  for  the  *P 
state,  204  space  and  306  space/spin  configurations  for 
the  * S  state,  and  42  space  and  42  space/spin  configura¬ 
tions  for  the  lS  state.11  Within  the  frozen  core  approxi¬ 
mation  the  accuracy  of  the  results  obtained  from  the 
full  Cl  calculations  is  only  limited  by  the  completeness 
of  the  basis  set . 

C.  Hartree-Fock  calculations  on  GaKr 

The  reference  orbitals  lor  the  POL-C1  calculations18 
on  the  GaKr  molecule  were  obtained  from  an  HF  calcu¬ 
lation  on  the  1  *2*  state  with  configuration11 

•••  13a*14a*15o*16a7v4  .  (la) 

In  the  POL-CI  calculations  the  core  orbitals  (12  o  orbit¬ 
als,  12  ir  orbitals,  and  4  6  orbitals)  are  always  required 
to  be  doubly  occupied  and  so  it  is  convenient  to  renum¬ 
ber  the  valence  orbitals  so  that  (la)  becomes 

1*2*:  la  *20*30  *401**  .  (lb) 

As  R  -  *>  the  above  orbitals  become 

la-4*,,  U~4ptKr 

2o~4poKr 

3a -4s0« 

4a  -  4paQ. 

The  1VO  method,*  with  a  la*2o*3o*l*4  core,  was  used 
to  generate  the  valence  It  and  Rydberg  5a  orbitals.  As 
R-»  these  orbitals  become  the  Apt  and  5s  orbitals  of 
the  gallium  atom.  The  virtual  orbitals  were  also  ob¬ 
tained  from  the  I  VO  calculations. 

In  terms  of  the  orbitals  defined  in  this  way  the  refer¬ 
ence  configurations  for  the  other  states  of  interest  are1* 

l*n:  la*2a*3o*l*42v  (2) 

2*2*:  la*2a*3a*5al»4  (3) 

llr:  la*2a*3a*U4  .  (4) 

D.  Polarization  configuration  interaction 
calculations  on  GaKr 

The  POL-CI  calculations1*  on  GaKr  Included  all  single 
and  double  excitations  relative  to  the  reference  config¬ 
urations  given  above  with  the  restrictions  that 

(1)  all  core  orbitals  remain  doubly  occupied  and 

(2a)  no  more  than  one  electron  occupy  the  Rydberg  5a 
end  virtual  orbitals  (*fl  and  *2*  states)  or 

(2b)  no  more  than  one  electron  occupy  the  Rytberg 
So  orbital  and  no  more  than  one  electron  occupy  the 
virtual  orbitals  (*Z*  states). 


For  the  *2*  states  the  calculations  considered  both  states 
simultaneously.  This  procedure  results  in  556  space 
and  1565  space/spin  configurations  for  the  *fl  states, 

764  space  and  2314  space/spin  configurations  for  the 
*2*  states,  and  368  space  and  556  space/spin  configura¬ 
tions  for.  the  '2*  state . 

The  less  restrictive  condition,  (2b),  is  necessary  for 
the  *2*  states  to  allow  for  configurations  such  as 

la*2a*[2wiv]5alir4  ,’ 

which  are  important  in  the  description  of  the  Rydberg 
*2*  state.-  As  was  the  analogous  case  in  the  gallium 
atom,  these  configurations  are  necessary  to  allow  for 
the  inclusion  of  the  4s*-4/>*  near-degeneracy  effect  in 
the  molecular  wavefunction. 

E.  Inclusion  of  spin-orbit  coupling  in  GaKr 

As  in  our  earlier  calculations  on  excimer  systems ,*’ 13 
we  have  adopted  a  simple  model3  for  including  the  effects 
of  spin-orbit  coupling  on  the  calculated  potential  energy 
curves  and  wavefunctions .  The  experimental  spin-orbit 
parameters  for  the  open-shell  atom,  gallium  in  the 
present  case,  are  used  to  determine  the  matrix  ele¬ 
ments  of  the  spin-orbit  Hamiltonian,  H„,  coupling  the 
molecular  states,  at  R  =•» .  These  matrix  elements  are 
then  assumed  to  be  independent  of  R  and  are  added  to 
the  diagonal  matrix  of  the  electronic  energies 

H°(fi)  =  6M«<(R)  +  H.  .  (5) 

The  energies  and  wavefunctions  with  spin-orbit  correc¬ 
tions  are  obtained  by  diagonalizing  H°(/t).  This  proce¬ 
dure  is  reasonable  only  so  long  as  (1)  the  molecular 
states  retain  the  identity  of  the  atomic  states  from  which 
they  arise  and  (2)  the  atomic  contributions  to  the  molecu¬ 
lar  spin -orbit  Interactions  are  dominant. 

For  GaKr  the  spin-orbit  interaction  affects  only  those 
Btates  which  arise  from  the  Ga(*P)  +  Kr(JS)  limit.  The 
Hamiltonian  matrix  for  the  ft  =  1/2  states  arising  from 
this  limit  is 


(6a) 

and  for  the  0  =  3/2  state  is 

H‘M'*«E(l*n)  +  X<J.  . 

(6b) 

In  (6)  \j,  is  one  third  of  the  *P, /*-*/* ./t  splitting  in  the 
gallium  atom.  The  energies  of  the  2*2'  and  1  *2*  states 
of  GaKr  are  unaffected  by  the  inclusion  of  spin-orbit 
coupling.  The  wavefunctions  (or  the  states  obtained  by 
diagonalizing  (6)  may  be  written  in  the  form 

|  / 1/0  «  cost)  1*11.0  ♦  ain&|  1  *2'o> 

(7a) 

|m/2>  -  -  sinfl|  1  *n.0  ♦  coe*|  1  *2 '«> 

(7b) 

|/S/2).|l*n.«)  , 

(7c) 

where  #  is  the  spin  rotation  angle.  The  wavefunctions 
hr  the  Mil  1/2)  and  I/© states  are 


(Wl/2)«|2*S*a) 
|/0-| I'E*)  , 


(7d) 

(7e) 
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TABLE  1.  Calculations:  on  the  2P,  ?S,  and  'S  (ion)  status  of 
Ga  and  the 's  state  of  krypton.  For  Ga  the  cx|>erimcntal  re¬ 
sults  have  been  corrected  for  spin-orbit  effects  (sec  the  text). 
Units  arc  as  Indicated. 


Atom  state  *P 

Ga(Ga’) 

‘s 

Kr 

*S 

Total  energies  (hartree)* 

POL-CI  0.22523 

0. 12130 

0.01912 

0.97377 * 

Full  Cl  0.23484 

0.12750 

0.02305 

... 

Excitation  energies  (eV) 

POL-CI  0.00 

2.  83 

5.61 

... 

Full  Cl  0.00 

2.92 

5.76 

... 

Expt'l'  0.000 

3.005 

5.930 

.  .  . 

Transition  moment  (eo0) 

POL-CI 

-  1 . 298<l 

... 

... 

Lifetimes  (nsec) 

POL-CI 

.  8.2* 

... 

... 

Expt’l 

6.8' 

•  •  • 

*  *  * 

'For  Ga  the  energies  are  relative  to  —  1923  hartree;  for  Kr  the 
energies  aro  relative  to  -  2751  hartree. 

bFor  Kr  the  POL-CI  wavefunction  is  equivalent  to  the  HF  wave- 
function. 

'Reference  4. 

*This  is  the  matrix  element  l/3l(*S|x  |,i>1)+  (*Sly  l*P,) 

»  <JS|rlJP.)|. 

•Using  the  experimental  excitation  energy  we  obtain  6. 9  nsec. 
'Reference  14. 


With  the  definitions  (7)  the  dipole  transition  moments 
coupling  the  III  1/2  state  with  all  of  the  lower  states  are 


p,(//f  1/2  - 1 1/2)  =  sinfl  (2*£*|  z|  1  *Z*)  (8a) 

M,(/// 1  /2  -  / 1/2)  =  cosfl <2 *£* |  x|  1  ’n^/ V5  (8b) 

p,(//7 1/2  -  7/ 1/2)  =  -  cosfl  <2  *E*|  z  1 1  *r)  (9a) 

ji,(//7  1/2  -  II 1/2)  =  sinfl  <2  *Z*|  x|  1  *n ,)//?  (flb) 

p,(OTl/2-/3/2)=-(2,z:*|x|l*n^//5  .  (10) 


Since  the  (2*E*lxll,n,)  transition  moment  is  expected 
to  be  comparable  in  magnitude  to  the  <2*z:*l*ll*E‘)  mo¬ 
ment  (at  R=«°  they  are  identical),  it  is  clear  from  (8) 
and  (9)  that  the  transitions  from  the  III  1/2  state  of  the 
/1/2  and  7/1/2  states  can  have  both  large  parallel  and 
perpendicular  components. 

IV.  RESULTS  FOR  Ga,  GaKr,  AND  GaKr' 

A.  Electronic  states  of  Ga,  without  spin-orbit 
corrections 

The  results  of  the  calculations  on  the  Ga  atom  are 
summarized  in  Table  I .  The  computed  *S—*P  excitation 
energy  ia  2.83  eV  (POL-CI)  and  2.92  eV  (full  Cl).  Av¬ 
eraging  the  multiplet  energies  for  the  *Pl/t  and  *P%n 
states  of  gallium  from  Moore,*  the  experimental  ‘S-'P 
splitting  la  calculated  to  be  3.005  eV,  Just  0.08-0.17 
eV  larger  than  the  computed  spacing-  The  errors  in  the 
calculated  ionisation  potentials,  5.61  sV  (POL-CI)  sad 
5.76  eV  (fciU  Cl),  are  somewhat  larger,  being  0.32  and 
0.17eV,  respectively. 


The  lifetime  of  the  *5  state  of  gallium  has  been  deter¬ 
mined  by  Norton  and  Gallagher14  to  be  6. 8 ±0.3  nsec. 

For  the  model  used  here  the  lifetime  of  the  *S  state  is 
independent  of  the  spin-orbit  corrections.  From  the 
POL-CI  wavefunctions  we  calculate  a  lifetime  of  8.2 
nsec  for  the  *S  state  (6. 9  nsec,  if  the  experimental  ex¬ 
citation  energy  is  used  instead  of  the  calculated  excita¬ 
tion  energy). 

B.  Electronic  states  of  GaKr  and  GaKr* ,  without 
spin-orbit  corrections 

The  energies  of  the  l*n,  1,2*2*,  and  1*2*  states  of 
GaKr  and  GaKr*  obtained  from  the  POL-CI  calculations 
are  listed  in  Table  II  and  the  resulting  potential  energy 
curves  are  plotted  in  Fig.  2.  As  is  usual  in  such  cal¬ 
culations,3**’ 13  the  energies  of  the  1*11  and  1*2*  states, 
both  of  which  arise  from  the  *P  limit,  are  not  exactly 
equal  at  R  =  15.0a„  (the  largest  value  of  R  considered). 
The  difference,  0.00114  hartree  (0.031  eV),  is  attribut¬ 
able  to  the  inequivalence  of  the  4/>o  and  4/>tt  orbitals  and 
to  core  polarization  effects  (the  core  orbitals  were  ob¬ 
tained  from  HF  calculations  on  the  1*2*  state  which  does 
not  have  the  full  rotational  symmetry  of  the  atom).  In 
the  plots  the  asymptotic  energies  of  the  states  have  been 
adjusted  to  give  the  experimental  atomic  energy  split¬ 
tings. 

In  line  with  the  discussion  in  Sec.  II,  the  potential  en¬ 
ergy  curve  for  the  1*11  state  is  found  to  be  less  repulsive 
than  that  of  the  1*2*  state,  thus  making  the  l!n  state  the 
ground  state  of  the  system.  In  fact,  we  find  that  the 
I*n  curve  is  slightly  bound,  £>#~0.04  eV  (see  Table  in). 
Although  spurious  minima  have  been  found  in  previous 
calculations  on  excimer  systems*’ 13  and  attributed  to 
basis  set  limitations,  the  well  in  the  1*D  curve  is  sub¬ 
stantially  larger  than  has  been  observed  heretofore.  We 
thuB  suspect  that  the  minimum  in  the  1  *11  curve  is  not 
just  a  result  of  calculational  limitations.  The  depth  of 
the  well  in  the  real  l*n  curve  is,  of  course,  expected 
to  be  significantly  larger  than  that  calculated  here  since 
the  POL-CI  method  is  not  designed  to  account  for  the  at¬ 
tractive  van  der  Waals’  interaction. 

As  predicted  in  Sec.  n  both  the  2*2’ and  1  *2*  curves 


TABLE  II.  Energies  obtained  from  the  POL-CI  calculations 
on  the  tow-lying  electronic  states  of  GaKr  and  GaKr*.  Dis¬ 
tance!  arc  in  bohr;  energies  are  in  hartree.  Energies  are 
relative  to  -  4674  hartree . 


R 

i*n 

GaKr 

1*E* 

2*£* 

GaKr* 

I  ‘r* 

3.75 

-1.02539 

-0.96866 

-  0. 93379 

-0.85625 

4.00 

-1.07933 

-1.01870 

-0.99238 

-  0. 90683 

4.50 

-1.14010 

-1.08688 

-1.06440 

-  0. 96758 

5.00 

- 1 . 17924 

-1.13516 

-1.08855 

-0.99287 

5.50 

-1.19225 

-1.16324 

-1.09545 

-1.00070 

6.00 

-1.19772 

-1.17938 

-  1 . 09677 

-  1.00202 

6. SO 

-1.19984 

-1.18665 

-  1 . 09631 

-  1.00107 

7.00 

-1.20045 

-1.19396 

-  1 . 09545 

-  0. 99950 

8.00 

-1.20020 

-1.19865 

- 1 . 09420 

-0.  99670 

10.00 

-1.10033 

-1.20029 

-1.09406 

-0.99402 

15.00 

-1.19970 

-1.20012 

- 1 . 09502 

-0.  99304 
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THE  LOW-LYING  STATES  OF  GoKr  AND 
GaKr* 


R(A) 


FIG.  2.  Calculated  potential  energy  curves  for  the  states  of 
GaKr  and  GaKr*  arising  from  the  Ga(,/>,  1S)  t  Kr('S)  and  Ga*('S) 

+  Kr(‘S)  separated  atom  limits.  The  curves  have  been  uniform¬ 
ly  shifted  to  correct  for  the  errors  in  the  gallium  atom  excita¬ 
tion  energies. 

are  found  to  be  bound  with  calculated  dissociation  ener¬ 
gies  of  0.047  and  0.24  eV,  respectively.  Again,  inclu¬ 
sion  of  the  attractive  van  der  Waals’  interaction  would 
be  expected  to  significantly  increase  the  calculated  well 
depths .  There  is  a  slight  hump  in  the  calculated  poten¬ 
tial  curve  of  the  2*2*  state,  -0.026  eV,  since  the  long- 
range  interaction  of  the  excited  Ga  and  Kr  atoms  is  re¬ 
pulsive.  This  hump  could,  however,  disappear  when  the 
van  der  Waals’  attraction  is  added  to  the  calculated 
curve. 

The  dipole  transition  moments  radiatively  coupling 
the  states  of  GaKr  are  given  in.  Table  IV  and  plotted  in 
Fig.  3.  At  R  =  “  the  transition  moment  coupling  the  1  *n 
and  1  *2*  states  is  Identically  zero  and  those  coupling 
the  (l*n,  1*2*)  states  with  the  22s*  state  are  approxi¬ 
mately  equal  to  the  *S-2JP  atomic  transition  moment 
(after  properly  accounting  for  differences  in  degeneracy 


TABLE  IB.  Spectroscopic  constants  for 
the  bound  ststes  of  and 

“oa^Kr*.  Units  are  as  indicated. 


i*n 

GaKr 

2lr* 

GaKr* 

1  V 

T„  eV 

0.00 

2.82 

3.40 

k 

3.73 

3.17 

3.14 

D,,  eV 

0.041 

0.047 

0.24 

u.*ft  cm*1 

36 

66 

83 

B.  0.0312  0.0442  0.04S2 


TABLE  IV.  Dipole  transition  moments 
coupling  the  low-lytng  stales  of  GaKr  ob¬ 
tained  from  the  POL-Cl  calculations. 


R 

iV-i2n*  2V-iiii*  2V-1V 

3.75 

-0.0037 

-0.7877 

0.4095 

4.00 

-0.8462 

-0.8680 

-0.3055 

4. so 

-0.1365 

-1.2090 

-2.0676 

s.oo 

0.0992 

-  1.2349 

-  1.4402 

5.50 

0.1257 

-1.2469 

-1.2456 

6.00  . 

0.1166 

-  1.2554 

-1.1893 

6.50 

0.0984 

-  1.2619 

-1.1800 

7.00 

0.0789 

-1.2677 

-1.1902 

8.00 

0.0476 

-1.2777 

-  1.2310 

10.00 

0.0163 

-  1.2898 

-1.2976 

15.00 

0.0018 

-1.2917 

-1.3101 

‘The  matrix  element  given  is  (n  V  lx  tl  *!!,>. 


factors).  As  R  decreases  rather  minor  changes  occur 
until  R~  2.5  A.  For  R<  2.5  A  substantial  changes  are 
noted  in  all  of  the  transition  moments,  although  the 
change  in  the  222*-l2n  moment  is  less  dramatic  than 
for  the  222*-l22*and  l22*-l2n  moments.  The  erratic 
behavior  of  the  transition  moments  for  R <  2.5  A  is  one 
manifestation  of  the  strong  interaction  of  the  1 22*  and 
2*2*  states  at  short  R. 

Although,  as  noted  above,  substantial  changes  are 
found  in  the  transition  moments  for  R <  2.5  A,  such  be¬ 
havior  can  be  expected  to  have  little  effect  on  the  ob¬ 
servable  properties  of  the  system .  At  R  =  2. 5  A  the  en¬ 
ergy  of  the  2*2*  state  is  >0.5  eV  above  its  asymptote 
so  that  the  region  R<  2.5  A  would  be  thermally  inacessi- 
ble. 


DIPOLE  TRANSITION  MOMENTS  AMONG 
THE  LOW-LYING  STATES  OF  GoKr 


R(A) 


FIG.  3.  Calculated  dipole  transition  moments  for  the  1  V- 
J*n,  2,E*-l,n>  and  2!S*-lV transitions  in  GaKr. 
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THE  LOW 'LYING  STATES  OF  GaKr  ANO 
GaKr*  WITH  SPIN'ORBIT  CORRECTIONS 


R(A) 

FIG.  4.  Calculated  potential  energy  curves  for  the  states  of 
GaKr  and  GaKr'  arising  from  tho  Ga(!jPty2,  j/2>  2Si/j)  +  Kr('Sj) 
and  Ga*('S0)*  K r ( !S0)  separated  atom  limits.  The  curves  have 
been  uniformly  shifted  to  correct  for  the  errors  in  the  gallium 
excitation  energies. 

C.  Electronic  states  of  GaKr  and  GaKr* ,  with  spin- 
orbit  corrections 

The  energies  of  the  states  of  GaKr  and  GaKr*  with 
6pin-orbit  corrections  are  given  in  Table  V  and  the  re¬ 
sulting  potential  energy  curves  are  plotted  in  Fig.  4. 

In  calculating  the  spin-orbit  corrected  energies  we 
shifted  the  1 2fl  energies  to  agree  with  the  1  lE*  energy 
at  the  largest  value  of  R  considered  (R  =  15.0aa).  Again, 
in  plotting  the  potential  energy  curves,  the  asymptotic 
energies  of  the  states  have  been  adjusted  to  give  the  cor¬ 
rect  atomic  energy  splittings. 

As  predicted  in  Sec.  IX,  at  short  R  the  curves  for  the 

TABLE  V.  Calculated  energies  of  the  1  1/2, 

11  l/2,  and  I  3/2  states  of  GaKr  with  spin- 
orbit  corrections.  Distances  are  in  bohr; 
energies  are  in  hartree.  Energies  are 
relative  to  the  energy  of  the  1 J2*  state  at 
R  -15  a,. ' 


R 

/  1/2 

/  3/2 

11  1/2 

3.75 

0.17228 

0.17484 

0.23152 

4.00 

0.11835 

0.12091 

0.18147 

4.50 

0. 04957 

0.05213 

0.11329 

5.00 

0.01841 

0.02099 

0.06502 

5.50 

0.00538 

0.00799 

0. 03698 

6.00 

-0.00015 

0.00252 

0. 02089 

6.50 

-0.00234 

0. 00040 

0.01170 

7.00 

-0.00307 

-0.00022 

0. 00650 

8.00 

-0.00316 

0.00003 

0.00215 

10.00 

-  0. 00280 

0.00090 

0.00103 

15.00 

-0.00251 

0.00126 

0.00126 

TABLE  VI.  Calculated  dipole  transition  moments  coupling  the 
111  1/2  and  1  1/2,  II  1/2  and  I  3/2  states  of  GaKr  with  spin- 
orbit  corrections.  Distances  arc  in  bohr,  moments  arc  in 
atomic  units. 


III  1/2- 

-/  1/2 

111  1/2-11  1/2 

III  1/2 -I  3/2 

R  M 

(x.y) 

z  (x.y) 

U.y) 

4.00  -  0.0086 

0.6135 

-0.3054  -0.0173 

-0.6138 

4.50  -0.0576 

0.  8545 

-2.0667  -0.0238 

-  0.8549 

5.00  -0.0549 

0.  8726 

-  1.4393  -0.0333 

-0.8732 

5.50  -0.0701 

0. 8803 

-  1.2437  -0.0496 

-0.8817 

6.00  -0.1007 

0.8845 

-  1.1850  -0.0752 

-0.6877 

6.50  -0.1504 

0.  8850 

-  1.1704  -0.1138 

-0.8923 

7.00  -0.2249 

0.  8803 

-  1.1688  -0.1694 

-0.8964 

8.00  -0.4405 

0.  8436 

-  1.1495  -0.3233 

-0.9035 

10.00  -0.7254 

0.7562 

-1.0759  -0.5098 

-0.9120 

16.00  -0.7564 

0.7457 

-  1.0697  -0.5273 

-0.  9134 

71/2  and  73/2  states  become  nearly  degenerate  and  just 
represent  the  curves  for  the  two  spin-orbit  components 
of  the  1 2n  state.  Also,  at  short  R  the  771/2  curve  close¬ 
ly  approximates  the  curve  for  the  12e*  state.  The  po¬ 
tential  energy  curves  for  the  22E*  (7771/2)  and  1  '£*(/ 0) 
states  are  unchanged  by  the  inclusion  of  spin-orbit  cor¬ 
rections. 

The  dipole  transition  moments  radiatively  coupling  the 
states  of  GaKr  with  spin-orbit  corrections  are  given  in 
Table  VI  and  are  plotted  in  Fig.  5.  It  should  be  noted 
that  both  the  z  component  of  the  7/7 1/2-7 1/2  transition 
and  the  x  component  of  the  7771/2-7/1/2  transition  are 
now  found  to  vary  significantly  even  for  R  >  2. 5  A.  This 
is  due  to  the  changing  nature  of  the  71/2  and  7/1/2  states 


DIPOLE  TRANSITION  MOMENTS  AMONG 
THE  LOW-LYING  STATES  OF  GaKr  WITH 
SPIN- ORBIT  CORRECTIONS 


R  ( A  > 

FIG.  5.  Calculated  dipole  transition  moments  for  the  III  1/2- 


71/2,  777  1/2-7  3/2,  and  III  1/2-7/  1/2  transitlona  in  GaKr. 
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TABLE  VU.  Excitation  energies  and  ion 
n  ization  potentials  for  the  Ga,  In,  and  T! 
atoms,  in  eV.  Taken  from  Itcf.  4. 


State 

Ga 

Group  IIIA  atoms 

in 

Tl 

’*1/1 

0.000 

0.  000 

0.000 

'P>n 

0. 102 

0.274 

0.966 

*sut 

3.073 

3.022 

3.282 

lc 

*0 

5.998 

5.786 

6.108 

as  the  atomic  spin-orbit  coupling  is  quenched  by  molec¬ 
ular  formation. 

V.  EXTRAPOLATION  TO  InKr  AND  TIKr 

To  provide  information  on  the  InKr  and  TIKr  mole¬ 
cules,  the  latter  being  the  most  experimentally  accessi¬ 
ble  of  the  Group  IIIA -krypton  molecules,  the  potential 
curves  for  InKr  and  TIKr  have  been  estimated  from  the 
GaKr  curves .  The  excitation  energies  and  ionization 
potentials  for  the  series  Ga,  In,  and  TI  are  given  in  Ta¬ 
ble  VO.  As  can  be  seen,  this  series  does  not  form  a 
steady  progression:  In  has  a  lower  ionization  potential 
and  excitation  energies  than  Ga,  as  expected,  but  Tl  has 
a  higher  ionization  potential  and  excitation  energies . 
This  anomalous  behavior  in  Tl  is  due  in  part  to  the  fill¬ 
ing  of  the  4/  shell  f' lanthanide  contraction")  and  to  the 
larger  spin-orbit  interactions.*5  These  effects  are  only 
partially  accounted  for  in  the  present  models. 

To  simulate  InKr  and  TIKr,  the  experimental  spin- 
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FIG.  7.  Model  potential  energy  curves  for  the  states  of  TIKr 
and  TIKr*  arising  from  the  Tl(JP,/2t  J/2,  *  Krl'i’,,)  and 

ln*(‘S0)  +  Kr(*S0)  separated  atom  limits,  see  the  text. 

orbit  parameters  for  In  and  Tl,  see  Table  VII,  were 
used  to  couple  the  curves  calculated  for  GaKr.  The  po¬ 
tential  energy  curves  for  InKr  and  TIKr  obtained  in  this 
tyay  are  expected  to  be  qualitatively  correct;  thp  curves 
are  plotted  in  Figs.  6  and  7.  As  before,  the  plotted 
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FIG.  6.  Model  potential  energy  curves  for  the  states  of  InKr 
and  InKr*  arising  from  the  In(JPi/j,  j/j,  *S|/j)  *  Kr('St)  and 
In*(*S,)»  Kr(*S,)  separated  atom  limits,  ace  the  text. 


DIPOLE  TRANSITION  MOMENTS  CONNECTING 
THE  LOW-LYING  STATES  OF  TIKr  WITH 
SPIN-ORBIT  CORRECTIONS 


1.5  2.0  2,5  3  0  3.5  4.0  4.5  SO  5  5  6.0 

RCA) 


FIG.  8.  Model  dipole  transition  moments  for  the  III  1/2-/  1/2. 
Ill  1/2-/  3/2,  and  III  1/2—//  t/2  transitions  in  TIKr. 
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FIG.  9.  Calculated  absorption,  AC„(T),  and  stimulated  emission. 


curves  have  been  shifted  to  give  the  correct  atomic 
excitation  energies.  The  effect  of  the  increasing 
spin-orbit  interaction  in  the  sequence  Ga<In<  T1  is  evi¬ 
dent  in  Figs.  4,  6,  and  7. 

Using  the  transition  moments  obtained  from  the  caicu  - 
lations  on  GaKr  and  the  wavefunctions  obtained  from  the 
TIKr  simulation,  the  transition  moments  of  TlKr  have 
been  estimated.  The  moments  so  obtained  are  plotted 
in  Fig.  0. 

VI.  ABSORPTION  AND  STIMULATED  EMISSION 
COEFFICIENTS  FOR  GaKr,  InKr,  AND  TIKr 

The  interest  in  the  Group  IDA-rare  gas  systems  arises 
from  the  possibility  of  their  use  as  visible  User  sys¬ 
tems.  In  order  to  judge  their  usefulness  as  lasers  it  is 
convenient  to  calcuUte  the  pure  absorption,  KV(T),  and 
stimulated  emission,  g,(T),  coefficients  for  the  perturbed 
atomic  transitions.  Obtaining  quantum  mechanical  re¬ 
sults  for  these  quantities  would  require  a  complex  cal¬ 
culation  which  is  not  justified  by  the  extrapoUtions  used 
to  obtain  the  InKr  and  TIKr  curves.  Consequently,  we 
have  used  instead  the  method  of  GalUgher  and  co- 
workers,1*  which  is  based  on  the  classical  Franck - 
Condon  principle.  In  this  approximation 

g,{T\,  =  [m*  lxXx78*K<  ^X^Ao)C(vT) 

K„(T)  =  (Mlxl(  f/g)g,(T)j  exp(h(v  -  »t)/kT) 
and 


1  V{RM}g„ 
di'/dR{v)  D(°°)  gA 


xexp{-  /[R(i/)]AT}  . 


In  these  equations,  J  refers  to  the  bands  associated  with 
the  transition,  uB  and  x0  are  the  frequency  and 

wavelength  of  the  atomic  transition  with  a  transition  rate 
of  <40(«I),  g*  *  2  for  the  *S1/Z  state  and  g-  2J  +  1  for  the 
*P j  states,  D{R)  is  the  transition  dipole  moment  at  R , 
gm  and  gA  are  the  statistical  weights  of  the  excited  mo¬ 
lecular  and  parent  atomic  state,  and  V*(R)  is  the  excited 
state  potential  curve  relative  to  the  energy  of  the  ex¬ 
cited  atomic  sUte.  [M],  [M*  ],  and  [x)  are  the  concen¬ 
trations  of  ground  and  excited  metal  atoms  and  of  rare 
gas  atoms  at  the  temperature  T. 


In  order  to  obtain  g„(T)  and  K,(T),  the  calcuUted 
curves  were  first  fit  with  cubic  splines .  The  spline  fit 
was  then  used  to  calculate  dv/dR,  and  these  quantities, 
along  with  the  atomic  transition  rates,1*  were  used  to  cal¬ 
culate  the  absorption  and  stimulated  emission  coeffi¬ 
cients  for  pressure  and  excitation  conditions  relevant  to 
the  experimental  studies.  We  have  calculated  the  coef¬ 
ficients  for  two  different  types  of  conditions.  The  high 
temperature  results  correspond  to  the  case  where  the 
concentration  of  the  metal  is  obtained  from  the  vapor 
pressure  of  the  metal  itself,  while  the  low  temperature 
results  correspond  to  obtaining  the  required  concentra¬ 
tion  of  the  metal  from  vaporization  of  MI,.  'This  latter 
condition  has  been  suggested  by  Gallagher2  as  a  possible 
means  of  obtaining  high  concentrations  of  the  metal  at 
low  temperatures.  In  both  cases  the  densities  used  are 
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PIG.  10.  Model  absorption,  K„(T),  and  stimulated  emission,  coefficients  for  InKr  at  r  =  300  and  1200  *K. 

1020/cms  for  [Kr|,  lO'Vcm5  =  3[M(,Pl/1)]  =  especially  accurate,  we  have  not  calculated  g„  and  K. 

and  2x  10‘Vcm1  /or  (M(*S1/2)].  Since  the  blue  wings  of  in  these  regions. 

the  lSt/z~*Pt/t,t/i  bands  are  due  to  transitions  occurring  The  resuitjng  absorption  and  stimulated  emission  co¬ 
at  large  ft,  where  our  curves  are  not  expected  to  be  efficients  for  GaKr,  InKr,  and  TIKr  are  given  in  Figs. 
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9-11.  In  these  figures  the  solid  lines  are  the  total  coef¬ 
ficients  calculated  by  assuming  that  the  transition  mo¬ 
ment,  D(R),  is  a  constant.  The  X’s  are  the  coefficients 
obtained  by  using  the  spatially  averaged  values  of  D(R) 
calculated  from  Table  VI  and  the  appropriate  spin-orbit 
parameters.  For  GaKr  at  1500 °K  (Fig.  10)  and  for 
403. 4<  X<  417.3  nm,  K„  >  10"3  cm'1.  The  dip  in  g„  close 
to  the  line  center  is  due  to  the  positive  value  of  V*(R 
~  6aa)-V* (R  =  «>). 

The  gain  coefficient,  which  is  approximately  equal  to 
g„-K,  can  be  estimated  from  these  curves.  For  the  fre¬ 
quencies  at  which  gain  occurs,  the  use  of  D(R)  decreases 
K„,  and  the  gain  by  30-4<ft  for  GaKr  and  by  9-15% 
for  TIKr.  For  TlKr  the  errors  introduced  by  the  ex¬ 
trapolation  procedure  are  likely  to  be  larger  than  those 
produced  by  using  D(«>)  rather  than  D(R).  The  maximum 
gain  occurs  approximately  at  the  minimum  in  the  /// 1  /2 
curve.  Because  depends  exponentially  on  V* ,  changes 
in  the  well  depth  of  the  excited  state  would  have  a  sig¬ 
nificant  effect  on  g„ .  If  the  well  depth  were  larger  tor 
TIKr,  as  Gallagher  and  co-workers  predict,*  g„  and  the 
gain  would  be  larger.  Likewise,  since  the  excited  state 
in  TIXe  is  predicted*  to  be  more  bound  than  in  TIKr,  the 
gain  should  be  larger  for  TIXe  than  for  TIKr . 
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Abstract 


It  is  demonstrated  with  the  help  of  model  calculations 
that  a  laser  field  can  be  used  to  supress  the  background 
in  electron  scattering  cross  sections  and  to  pick  out  only 
those  parts  which  are  rapidly  varying  as  function  of  the 
incoming  electron  energy. There fore  it  becomes  possible  to 
measure  the  pure  Breit-Wigner  peaks  and  threshold  effects 
independent  of  the  fact  that  interference  between  the 
resonances  and  the  background  in  radiationless  electron 
scattering  is  strong. The  results  are  interpreted  within 
the  low  frequency  limit  for  free-free  transit ions. The 
AC  Stark  shift  of  the  resonance  is  also  observed. 
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1 .  Introduction 

Some  of  the  most  spectacular  physics  in  electron-atom 
scattering  is  in  the  appearance  of  resonances  and  threshold 
effects. Unfortunately, in  many  cases  these  effects  are  small 
structures  hurried  in  a  big  background  and  their  shape  is 
usually  modified  by  interference  with  the  background. 
Sometimes  it  is  difficult  to  be  sure  if  a  measured  structure 
is  real  or  noise  and  sometimes  it  is  hard  to  pull  the 
interesting  physical  parameters  out  of  the  measured  data. 

For  a  review  of  these  effects  and  their  experimental  diffi¬ 
culties  see  ref.1  .If  on  the  other  hand, an  experimentalist 
were  able, by  some  method, to  supress  all  the  background  in 
the  neighborhood  of  a  resonance, then  a  pure  Breit-Wigner 
peak  would  be  observed  and  it  would  be  straightforward  to 
determine  the  position  Eg  and  the  width  T  of  this  resonance. 
The  purpose  of  this  paper  is  to  explain  a  method  of  perfor¬ 
ming  an  electron  scattering  experiment  in  such  a  way  that 
only  those  parts  of  the  scattering  amplitude  lead  to  a 
signal, which  vary  rapidly  as  function  of  the  incoming 
electron  energy. Thereby  we  pick  resonances  and  cusps  out 
of  the  total  scattering  amplitude  and  get  rid  of  all 
background  influences. 

The  basic  idea  is  to  perform  the  electron  scattering 
process  within  a  strong  laser  field  and  to  collect  only 
those  scattered  electrons  whieh  have  emitted  or  absorbed 
a  certain  number  of  photons. This  method  has  been  suggested 
first  in  ref .2  in  a  treatment  of  free-free  transitions  in 
first  order  perturbation  theory  and  low  frequency  limit. 
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It  will  become  clear  in  this  paper  that  the  idea  works 
under  much  more  general  circumstances. Processes  in  which 
an  electron  scatters  off  an  atom  or  molecule  and  emits  or 
absorbs  photons  at  the  same  time  are  known  under  the  name 
free-free  transit ions. (For  a  review  of  these  processes 
see  ref.3)*In  these  processes  three  systems  (  electron, 
target, laser  field  )  interact  with  each  other  and  this 
makes  a  theoretical  treatment  so  complicated  that  up  to 
now  only  special  cases  and  approximations  have  been 
investigated. If  the  photon  energy  is  far  away  from  all 
transition  energies  of  the  target  states  involved  in  the 
process  then  the  interaction  between  the  laser  and  the 
target  can  to  a  high  degree  of  approximation  be  neglected 
and  the  problem  reduces  to  the  description  of  electrons 
under  the  simultaneous  influence  of  the  target  potential 
and  the  laser  field. To  check  this  the 

»—  influence  of  the  laser-target  interaction  on  the 

processes  considered  here  can  be  estimated  by  the  following 
considerations: The  external  field  induces  an  electric 
dipole  moment  in  the  target  atom  and  this  dipole  field 
can  be  felt  by  the  scattered  electron  and  can  thereby 
modify  the  electron-target  interaction  potential  V.Next 
let  us  estimate  the  order  of  magnitude  of  this  effect. 

In  most  of  our  calculations  a  photon  energy  of  0.003  a.u. 
is  used  which  is  close  to  the  COg  laser  photon  energy  of 
117s»V. Because  this  photon  energy  is  smaller  by  a  factor 
10  to  100  compared  to  the  first  excitation  energy  of  the 
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atomic  ground  states  we  use  the  static  polarizability  to 
calculate  the  induced  dipole  moment  (  Thereby  we  overesti¬ 
mate  the  effect  ).For  atoms  the  polarizability  is  between 
2*10”2'’  cm^  (  in  e.s.  units  )  for  such  a  small  closed 
shell  atom  like  Helium  and  400* 10”2^  cm^  for  such  a  big 
alkali  atom  like  Cesium. 

t-  To  see  strong  signals  in  free-free  transitions  we  need 
laser  power  fluxes  in  the  order  of  10  Watt/cm  and  for  a 
CO2  laser  this  corresponds  to  an  electric  field  amplitude 

of  270000  V/cm  or  900  e.s.u. .Therefore  the  induced  dipole 

-22  -20 

moment  is  between  1.8*10  e.s.u.  and  5.6*10  e.s.u. 

for  the  various  atoms. A  scattering  electron  at  a  distance 
—8 

of  10  cm, which  is  in  the  order  of  magnitude  of  the 

average  distance  of  an  electron  in  a  resonance  state  from 

the  atomic  center, will  feel  a  potential  energy  Eg  of  a 

_2  * 

magnitude  between  10  ^eV  and  10  eV.In  electron  scattering 
dipole  effects  are  usually  of  importance  only  at  very 
low  energy  where  the  electron-dipole  interaction  energy 
Eg  is  of  the  same  order  of  magnitude  as  the  kinetic  energy 
of  the  electron. In  the  case  considered  here  Eg/E^&O.OI 
in  the  worst  case .There fore  it  seems  to  be  justified  to 
neglect  the  laser-atom  interaction. 

This  neglect  of  the  laser-target  interaction  does  not 
of  course  apply  at  all  to  molecules  which  have  vibrational 
transition  energies  close  to  the  photon  energy. In  general, 
the  laser  frqquency  must  be  far  out  of  resonance  with  the 
target  system. 
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In  addition  the  laser  field  strength  considered  here 
is  small  compared  to  the  internal  field  in  the  atoms  and 
therefore  it  is  not  necessary  to  worry  about  the  ionization 
of  the  target. 

Now  even  this  simplified  problem  without  laser-target 
interaction  has  not  yet  been  solved  in  general. Various 
treatments  have  either  treated  the  laser  field  in  pertur¬ 
bation  expansion  or  have  used  the  low  frequency  approxi¬ 
mation. The  method  used  here  to  supress  the  background  can 
be  interpreted  theoretically  in  the  low  frequency  limit. 

The  relevant  ideas  of  this  approximation  will  be  given 
later  in  section  3  (  For  some  additional  information  about 
more  general  aspects  of  the  low  frequency  approximation 
see  refs.  4-7  ). 

To  fully  appreciate  the  method  it  is  perhaps  better  to 
forget  about  formal  theories  and  instead  to  make  some 
simple  model  calculations  in  the  spirit  of  a  numerical 
experiment. In  section  2  we  investigate  in  this  way  what 
may  happen  in  an  experiment  since  up  to  now  the  new  method 
has  not  yet  been  realized  in  the  laboratory. As  such  we 
construct  a  simple  1-dimensional  model , which  can  be  solved 
in  terms  of  known  functions. We  do  not  use  the  low  frequency, 
approximation  and  we  do  not  make  a  perturbation  expansion 
of  the  electron-laser  interaction. The  approximations  we 
must  make  are  the  neglect  of  the  laser-target  interaction, 
mentioned  above ,the  dipole  approximation  and  a  cut  off  of 
the  photon  Pock  space  due  to  the  inability  to  numerically 
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handle  more  than  a  finite  number  of  photon  number  states. 

In  practice  this  is  not  a  problem  as  all  reported  results, 
upon  inclusion  of  further  photon  states  in  the  model, do 
not  significantly  change. 

In  section  3  we  present  a  theoretical  interpretation 
of  our  numerical  results. Section  4  contains  conclusions 
and  final  remarks  and  some  suggestions  onto  how  to 
optimise  possible  experiments. In  the  appendix  we  explain 
our  model  in  detail  and  show  how  we  calculate  all  transition 
amplitudes  in  terms  of  known  functions. 

We  use  atomic  units  for  all  numbers  throughout  the 


whole  rest  of  the  paper 
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2.  Results  of  the  numerical  experiments 

We  take  a  scattering  problem  in  1  dimension  and  represent 
the  target  by  a  2-state  system  with  an  excitation  energy 
E12=0-6  .We  assume  that  this  target  provides  a  2-channel 
square  well  potential  V  to  the  incoming  electron  (  see  fig.1). 
In  each  figure  below  we  have  indicated  the  special  choice  of 
V  used. We  always  start  the  process  with  the  target  in  its 
lower  state. If  the  laser  is  switched  off  and  the  incoming 
kinetic  energy  E^  of  the  electron  is  below  0.6  then  the 
electron  may  either  be  transmitted  or  reflected  without 
energy  change. We  denote  this  probability  for  transmission 
by  T  .  and  the  probability  for  reflexion  by  R  .  (  rl  is  an 
abbreviation  for  "radiationless"  ).If  E^  is  above  0.6  the 
electron  may  excite  the  target  and  leave  it  in  its  upper 
state .We  denote  the  probability  for  excitation  in  radiation¬ 
less  transmission  by  Srl  and  the  probability  for  excitation 
in  radiationless  reflexion  by  Srl  and  Qrl  exist  only 

for  Ei>  0.6  . 

If  we  switch  on  the  laser  field, then  the  electron  can, 
in  addition  to  the  above , exchange  energy  with  the  field  in 
integer  multiples  of  the  photon  energy  u>.We  denote  the 
probabilities  for  processes  in  which  the  electron  gains 
( loses )  the  energy  Nu from  the  field  by  Qn,RN,ER’TN 
respectively. Rjj(E^)  is  the  probability  that  an  electron 
comes  in  with  initial  kinetic  energy  E^,hits  the  target, is 
reflected  leaving  the  target  behind  in  its  lower  state 
and  flies  away  with  a  final  kinetic  energy  E^ ■  E^  +  N  w  . 
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Similarly ,Qjj(E^)  is  the  probability  that  the  electron 
comes  in  with  energy  E^,hits  the  target, is  reflected  leaving 
the  target  behind  in  its  upper  state  and  flies  away  with 

final  energy  Ef=Ei  -Ei2  +NtJ*Tu^Ei^  and  SN^Ei^  are  tiie 
probabilities  for  the  corresponding  processes  in  transmission. 
Of  course  ,Qn(E^)  and  Sjj(E^)  only  exist  for  Ef>0  i.e. 

E^  +  Nco  >  E^2  -The  quantities  Q,R,S,T  are  the  1-dimensional 
analogs  to  the  differential  cross  sections  in  3  dimensions. 

In  fig .2  we  have  chosen  a  potential  V  which  causes  an 
elastic  scattering  resonance  at  an  incoming  energy  of 
Er^0.233  with  a  width  of  T W  0.0002. This  P  is  in  the 
same  order  as  the  natural  width  of  electron-rare  gas  atom 

is  far  below  the  threshold  at  ET=0.6 
and  therefore  no  excitation  of  the  target  is  possible  near 
this  energy. In  the  top  line  of  fig.2  we  show  the  quantities 
Trl(Ei)  on  the  left  and  Rrl(E^)  on  the  right  as  functions 
of  E^. Below  we  show  the  quantities  T ( E ^ )  and  R^E^)  as 
functions  of  for  several  values  of  N.The  photon  energy 
is  always  CO  *0.005  which  is  close  to  the  photon  energy  of 
a  CO^  laser. The  power  density  of  the  field  is  so  big  that 
the  quantity  o< «  eA  /  me  "fe w  has  a  magnitude  of  exactly  1. 

(  e  is  the  electron  charge, m  is  the  electron  mass,c  is  the 
speed  of  light, A  is  the  amplitude  of  the  vector  potential). 
This  is  still  a  moderate  power  density  and  therefore  only 
one  photon  processes  give  a  strong  signal. Two  photon 
transitions  can  Just  be  seen. In  fig. 3  we  show  the  results 
for  a  calculation  in  which  we  have  increased  the  laser 


resonances. The  energy  E^ 
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power  so  that  <x  =2  but  left  all  other  parameters  unchanged. 

The  important  results  of  these  figures  are: The  elastic 
resonance  causes  a  whole  series  of  resonances  in  each  R^ 
and  Tj.  but  only  a  few  ones  are  strong  enough  so  that  we  can 
observe  them  immediately .The  distance  between  two  adjacent 
structures  is  exactly  oj .The  strength  of  the  individual 
structures  depends  on  the  laser  power. In  general  we  see 
more  resonances  if  we  increase  the  laser  power. The  most 
striking  observation  is, that  in  TN,NfO  there  is  no  back¬ 
ground  at  all  and  the  resonances  appear  as  pure  Breit-Wigner 
peaks. In  all  other  cases  (  i.e.  for  TQ  and  all  )  the 
interference  of  the  resonances  with  the  background  is  either 
similar  to  the  corresponding  radiationless  case  only  with 
the  difference  that  the  relative  resonance  effect  is  smaller 
or  the  resonance  shape  is  reflected  compared  to  the 
radiationless  case. The  best  examples  for  this  reflection 
are  in  R+^  in  fig. 5. We  see  that  it  is  easier  to  determine 
Er  and  r  out  of  T+/)  or  T_^  than  out  of  Trl  or  Rr^*Tn  fig. 2 
note  that  Tjj  and  T_N  or  Rjj  and  R-N  look  similar  in  shape 
but  are  shifted  by  Nw  . 

In  fig. 4  we  have  changed  the  offdiagonal  elements  of  V 
so  that  the  width  T  of  the  elastic  resonance  is  now  0.0012 
which  is  in  the  order  of  the  resolution  of  common  electron 
spectrometers. This  change  in  V  causes  also  a  small  shift 
in  the  position  Ejj  of  the  resonance  (  see  the  top  line  in 
fig. 4  ).A11  other  parameters  are  the  same  as  in  fig. 2. The 
width  F  is  now  so  big  that  there  are  overlaps  of  the 


various  resonance  structures  in  Rjj  and  T^  which  influence 
the  resonance  shapes  in  general. For  T^,N+0  the  wings  of 
two  adjacent  resonance  peaks  simply  add  without  any  visible 
interference  .Again  and  T^  or  Rjj  and  R_^  look  similar 
in  shape  but  shifted  by  No). In  fig. 5  we  show  the  resonance 
curves  for  (X  =3  and  otherwise  the  same  parameters  as  in  fig. 4 
We  did  not  make  calculations  with  another  value  of  go 
because  the  interesting  quantity  is  the  ratio  between 
and  I”  and  a  increase  of  co  would  cause  similar  effects  as 
a  decrease  of  fT  . 

The  only  threshold  of  our  2-channel  model  is  at  E,p=0.6 
and  it  is  interesting  to  compare  radiationless  scattering 
and  free-free  transitions  close  to  this  energy. Fig. 6  shows 
radiationless  scattering  in  the  top  line  and  free-free 
transitions  below. In  the  drawings  T  and  R  are  represented 
by  solid  lines  and  S  and  Q  are  represented  by  dotted  lines. 
First  we  see  that  in  free-free  transitions  the  threshold 
for  Sjj  and  QN  lies  at  E^=E,j,  -  M  to.  This  is  easy  to  understand 
because  the  energy  which  the  electron  gains/losses  from 
the  laser  field  goes  into  the  energy  balance  for  the  excitati 
of  the  target  by  the  electron  impact  and  the  new  channel 
opens  as  soon  as  the  final  kinetic  energy  of  the  electron 
is  above  zero  i.e.  as  soon  as  E.  +  N  co  -  E^>  0  or 
E^  >  E^,  -  Nco.In  T_^  we  see  a  plateau  between  and  E^+co 
and  in  T+^  a  plateau  between  E^  -  CJ  and  E^,  with  a  sharp 
drop  on  both  sides. In  addition  there  are  smaller  threshold 
effects  shifted  by  co  away  from  the  sharp  drops. In  T_2 
and  T+2  we  see  a  sharp  peak  at  ET  +  co  or  E^  -  cO 


respectively  and  again  smaller  effects  shifted  away  by  w  . 
In  the  other  cases  (  TQ  and  all  RN  )  the  threshold  effects 
occur  at  various  +  nu>  and  their  shape  is  either  of 

the  same  type  as  in  the  corresponding  radiationless 
scattering  or  it  is  just  turned  upside  down. In  and 

Qn  there  are  only  very  weak  structures , too  weak  to  be 
seen  clearly  in  fig. 6. 

In  fig. 7  we  see  results  for  0(  =2  and  otherwise  the  same 
parameters  as  in  fig. 6. Now  the  threshold  effects  are 
distributed  over  more  energy  values  and  for  example  in 
T_^  the  main  drop  on  both  sides  of  the  plateau  is  shifted 
outwards  by  c*)  on  both  sides  compared  to  fig. 6. In  the 
picture  for  Qq  we  see  an  "accidential"  zero  near  E^=0.607 
which  will  be  explained  below. In  fig. 8  we  see  results 
for  u>  =0.002  and  otherwise  the  same  parameters  as  in  fig. 6. 
All  curves  look  nearly  the  same  as  in  fig. 6  only  the  energy 
scale  is  changed  by  a  factor  2.5  which  is  the  ratio  between 
the  photon  energies  in  the  two  examples. 

The  most  important  results  for  the  thresholds  are  in 
summary: A  sudden  drop  (  or  increase  )  in  the  radiationless 
scattering  as  function  of  the  incoming  energy  produces  in 
the  free-free  probabilities  T+y^  and  T_^  a  plateau  near  the 
threshold  energy  and  zero  signal  otherwise. At  moderate 
laser  power  the  length  of  the  main  structure  is  exactly  60  . 
Thereby  the  presense  or  absense  of  a  threshold  makes  a 
100%  change  in  the  relative  signal  strength  but  the 
absolute  signal  strength  is  always  very  low. Depending  on 
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the  values  of  oj  and  ot  the  strong  contributions  to  the 
threshold  effects  in  TN  appear  at  different  energies. 

3. Interpretation  of  the  numerical  results 
A  theoretical  explanation  of  the  results  of  section  2 
can  be  given  most  easily  within  the  low  frequency  appro¬ 
ximation  and  therefore  let  us  first  explain  its  basic  idea. 
In  order  to  be  more  general  we  give  all  formulas  for  the 
3-dimensional  case  and  consider  the  numerical  results 
from  section  2  as  the  special  cases  of  the  scattering 
angles  0  and  T T  .If  the  laser  wavelength  is  very  long 
compared  to  atomic  distances , then  we  can  separate  the  free- 
free  transition  into  three  steps. First  a  free  electron 
moves  within  the  laser  field  and  can  virtually  emit  and 
absorb  photons. Let  us  assume  that  the  laser  field  is 
single  mode  and  in  the  pure  number  state  J  in  absense 
of  the  electron. We  denote  the  state  of  a  free  electron 
moving  with  momentum  by  |  "j?  ^  .If  we  neglect  photon 
depletion  effects  and  use  the  dipole  approximation, then 
the  exact  state  of  the  electron  in  the  field  is  given  by 
(  see  appendix  ) 

!4>p*,N>  =  2.  Jn(c**P*  )l?>|N+n>  (1) 

J  is  the  Bessel  function  of  first  kind  and  order  n. 
n 

*  eA  ^  /  me  "feu*  where  ^  is  the  polarization  vector 
of  the  laser  field. This  is  the  3-dimensional  generali¬ 

zation  of  the  o(  given  above  in  section  2. 
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The  Bessel  functions  J  (  <*  p  )  can  be  viewed  as  interaction 
coefficients  for  virtual  n  photon  absorbt ion/emission  at 
absense  of  a  target. 

In  dipole  approximation  there  is  no  recoil  of  the 
electron  during  emission  or  absorbtion  of  a  photon  and 
therefore  the  electron  momentum  does  not  depend  on  how 
many  photons  the  electron  has  absorbed  or  emitted  i.e. 

|  p  y  in  (1)  is  independent  of  n.The  energy  of  the  electron 
is  different  in  its  various  states  and  the  electron  wave 
connected  with  the  photon  number  state  |  N+n^  is  at  energy 
E*p  /2m  -  ncj  .Only  the  n=0  term  is  on  the  energy  shell 
and  all  other  terms  are  off  shell. 

Then  in  the  second  step  of  the  free-free  transition 
this  mixture  of  electron  waves  hits  the  target  and  is  scattered. 
The  main  idea  of  the  low  frequency  approximation  is  to  neglect 
the  laser-electron  interaction  during  this  second  step 
and  therefore  the  electron-target  scattering  is  described 
by  a  radiationless  scattering  amplitude. But  according  to 
what  VOS  said  above  each  term  of  the  sum  in  (1)  is 
scattered  at  its  particular  intermediate  energy  E^  -  nco 
and  its  scattering  is  therefore  described  by  a  scattering 
amplitude  at  this  shifted  energy. 

In  the  third  step  the  scattered  electron  waves  interact 
again  with  the  laser  and  evolve  into  their  final  states. 

Along  these  lines  the  following  formula  (5)  has  been 
derived  first  in  ref. 8  and  confirmed  later  in  ref. 9  by 
another  derivation. This  formula  gives  the  scattering 
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amplitude  f^E^,-#’)  for  an  electron  to  come  in  with 
kinetic  energy  E^,be  scattered  by  an  angle  and  to  have  a 
final  kinetic  energy  E^  of 

Ef  =  Ei  +  N  co  (2) 

fN(Eit^)=r  JN_k(«  Pf)  frl(E.+ko,,^)  (3) 

f  1  is  the  corresponding  amplitude  for  radiationless 
electron-target  scattering. The  three  factors  correspond 
to  the  three  physical  stages  above. An  analoguous  formula 
holds  for  scattering  with  excitation  of  the  target  if  we 
take  on  the  r.h.s.  of  (3)  the  corresponding  scattering 
amplitude  for  radiationless  excitation  of  the  target  and 
take  the  excitation  energy  into  the  energy  balance  (2). 

Now  let  us  decompose  fy^  in  (3)  into  a  resonance  part 
and  a  background  part 

frl(Ei+k u  )-f®1(E.+k o  )+f®i(Ei+k  ~ .-S’ )  (4) 

The  Bessel  functions  decrease  rapidly  for  increasing  order 

and  fixed  argument  as  soon  as  the  absolute  value  of  the 

order  becomes  larger  than  the  absolute  value  of  the 

argument .Therefore  only  a  few  k  give  a  strong  contribution 

to  the  k  sum  in  (3)«  depends  only  weakly  on  its  energy 

BG 

argument  then  we  can  neglect  this  dependence , take  all  f^ 

•Dp 

at  the  energy  E^,pull  f^  out  of  the  k  sum  and  apply  the 
addition  theorem  of  the  Bessel  functions  and  find 

W^-V^Pf-Pi))  *??(*!.*>♦ 

lc 


Our  model  calculation's  independent  of  the  low  frequency 


approximation  and  so  a  comparison  between  the  exact 
numerical  results  of  sec. 2  and  (5)  can  be  viewed  as  test  of  (5)* 

The  present  state  of  the  experimental  verification  of  (5) 
is  as  follows:In  ref. 10  |fjj|^  bas  been  measured  as  function 
of  N  for  several  values  of  9  ,E^,"?.A11  results  are  in 
qualitative  agreement  with  (5). It  has  not  been  possible 
to  check  for  quantitative  agreement  because  the  exact  laser 
power  distribution  in  space  and  time  has  not  been  known. 

For  quantitative  calculations  its  knowledge  would  be 
absolutely  necessary  since  multiphoton  transitions  are 
a  nonlinear  effect  and  depend  essentially  on  the  exact 
power  distribution  in  space  and  time  and  not  only  on  the 
average  power  (  For  these  problems  see  also  ref. 11  ).In  ref. 10 
resonances  did  not  play  any  role  and  therefore  only  the 
first  term  in  (5)  has  been  measured. From  such  a  measurement 
we  cannot, in  contrast  to  the  present  work, learn  anything 
new  about  electron-atom  scattering  because  the  result  is 
essentially  the  product  of  the  elastic  background  cross 
section  and  a  Bessel  function  -  both  well  known  quantities. 

In  ref 8. 12, 13  jf^j^  bas  been  measured  as  function  of 
Ex  near  the  argon  resonances  at  11eV  but  unfortunately 
only  in  backward  diraction  for  the  one  scattering  angle 

160°. Since  a  weak  laser  has  been  used  only  the  resonance 
structures  at  and 
strength  is  in  full  agreement  with  (5). 

For  scattering  without  excitation  of  the  target  | 

since  <»)«,  and  therefore  in  foreward  direction  we  find 

■) 


Eg*  u>  could  be  seen. Their  shape  and 
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CX  (pf-p^)250.The  background  term  in  (5)  disappears  in 
foreward  direction  for  N$0  because  JN(0)=  Q.This  is  the 
explanation  why  we  did  not  see  a  background  contribution 

p 

for  Tjj,N$0  in  the  model  calculations. f  ^  produces  a 
resonance  structure  when  its  energy  argument  is  at 
and  therefore  we  get  a  resonance  structure  in  the  k-sum 
in  (5)  every  time  when  E^kw-Eg  for  a  k  a 2^ The  magnitude 
of  the  resonance  effect  at  E^=E^  -k  u>  depends  on  the  exact 
value  of  the  two  Bessel  functions  in  the  corresponding 
term  in  (5) -If  F<<co  we  get  a  series  of  well  separated 
resonances  with  an  energy  spacing  co  between  adjacent 
peaks. This  is  exactly  what  we  saw  in  figs. 2-5. 

In  backward  direction  p£*-p?  and  in  general  (p^-p^ 

is  not  near  a  zero  and  the  background  does  not  disappear. 
Depending  on  the  relative  signs  of  the  various  Bessel 
functions  the  relative  phase  between  the  resonance  term 
and  the  background  term  is  either  the  same  as  in  the 
radiationless  case  or  it  is  just  shifted  by  7 T  .Therefore 
the  resonance  shape  is  either  the  same  as  in  the  radiation¬ 
less  case  or  it  is  just  reflected. The  magnitude  of  the 
relative  resonance  effects  depends  again  on  the  exact 
values  of  the  corresponding  Bessel  functions  and 
thereby  on  the  parameters  E^CX.We  saw  that  TN  and  T_N 
or  fig  and  *-N  look  nearly  the  same  in  shape  but  are 
only  shifted  by  Nfcj.This  can  be  explained  by  the  fact 
that 

One  effect  which  is  not  contained  in  (5)  is  the  AC 
Stark  shift  of  the  resonances. Let  us  look  at  fig. 9. 
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In  this  figure  we  show  resonance  structures  in  Tg  and 

Rjj  at  various  initial  energies  E^=Eg+kco  and  for  various 

values  of  cu ,  o<  and  N.The  potential  V  is  the  same  as  in 

figs. 2  and  3*In  each  picture  we  show  the  low  frequency 

approximation  of  (3)  as  broken  line  and  the  exact 

numerical  result  as  solid  line. We  see  a  shift  of  the 

resonance  structures  which  does  not  depend  on  E^ ,N  and 

but  depends  on  cj  and  ot  .Besides  this  shift  the  low 

frequency  limit  gives  the  correct  shape  and  height  for 

all  resonance  peaks. The  fsct  that  the  shift  does  not 

depend  on  the  particular  peak  at  which  we  look  i.e.  does 

not  depend  on  E^,N  and  &  is  anindication  that  the 

resonance  itself  is  shifted  and  not  only  its  appearance 

in  the  free-free  transitions. The  Stark  shift  goes  quadra- 

tically  with  cj  and  quadratically  with  cx.This  is  in 

exact  agreement  with  the  results  of  ref. 9. In  our  examples 

P  P 

in  fig. 9  the  shift  AEg  «-0. 23  w  «  .According  to  ref. 9 

2  2 

the  proportionality  constant  between  AEg  and  u  «  is  a 
measure  of  the  laser  induced  coupling  between  the  resonance 
wave  function  in  channel  2  and  the  continuum  wave 
function  in  channel  1. Therefore  a  measurement  of 
AEfl/oj  oc  could  give  experimental  information  about 
the  dipole  matrix  element  ^<fl  I  l  • 

We  could  not  see  a  change  of  P  with  increasing  laser 
power. For  this  extra  width  caused  by  the  Stark  decay  to 
become  important  much  higher  laser  powers  would  be  necessary. 


Now  let  ua  try  to  explain  the  threshold  effects  in 
figs. 6-8  by  making  a  quite  crude  model  for  a  threshold 
effect  in  radiationless  scattering  and  representing  the 
amplitude  for  scattering  without  target  excitation  by 

frl(Ei’^)=c(^)+0(Ei“ET)  d(^)  (6) 

where  c  and  d  do  not  depend  on  and  @  is  the  unit 

step  function. If  we  insert  (6)  into  (5)  we  get 

fw(Ei,3’)-c(*9')  Jn(&(pJ-pJ))  + 

d(3  )  ^  JN_k(wp^)  ©(E^^+k  to -Et)  Jk(-}p*j)  (7) 

— ♦  — ^ 

For  p^j^p^  and  NfO  the  first  term  drops. The  second  term 
produces  a  step  at  each  E^=E^,-kco  ,k  e  ,2  .The  relative 
strength  of  the  various  steps  depends  on  the  exact  values 
of  the  Bessel  functions  and  thereby  on  the  values  of  the 
parameters  Ei  and  CX  • 

If  we  set  •&"  =0  in  (7)  and  assume  a  moderate  laser  power 
for  which  only  first  order  contributions  are  important, 
we  obtain 

f-^Ei»°>i  <*|pI  d(O)-[0(Ei-ET)-e(Ei-<O-ET)j' 

This  is  just  a  plateau  between  ET  and  Ey+  CJ  and  it  is  zero 
otherwise. An  analoguous  calculation  applies  to  T+/|  and 
gives  a  plateau  between  E^-co  and  Eij.T^  in  fig. 6  does 
not  have  a  sharp  step  at  E,j  but  a  rounded  one  and  therefore 
also  T+<1  and  T_^  have  a  rounded  drop  off  on  both  sides  of 
the  plateau. 

We  see  additional  smaller  threshold  effects  at  other 
Ej,  which  come  from  higher  order  contributions. For  higher 
laser  powers  as  e.g.  in  fig. 7  higher  order  contributions 
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may  become  more  important  than  first  order  contributions 
and  therefore  the  steps  at  other  E^  become  more  important. 

If  we  look  at  T_2  and  take  only  second  order  contri¬ 
butions  of  (7)  we  get 

f_2(Ei,0)*d(0)  5«p2£e(Ei-ET )-2d(Ei-ET- W)+  e(Ei-ET-2co)]> 
This  is  a  positive  plateau  between  E^  and  E^+  oo  and  a  negative 
plateau  of  the  same  absolute  value  between  E,p+u>  and 
Ej+2co  and  it  is  0  otherwise. T_2»pf  j  f_2(E^,0)j^  is  then 
just  a  plateau  between  E,p  and  ET+2«o.The  step  in  T  ^  in 
fig. 6  is  rounded  and  therefore  a  wedge  remains  at  E,j,+  cj 
in  T_2  in  fig. 6. This  wedge  is  a  left  over  between  two 
rounded  cut  offs  on  both  sides. For  higher  power  density  in 
fig. 7  again  additional  threshold  effects  at  other  E^  become 
important  in  T_2.An  analoguous  reasoning  applies  to  T+2. 

For  Tq  and  both  terms  in  (7)  contribute  and  the 
relative  phase  between  both  terras  is  either  the  same  as 
in  the  radiationless  scattering  or  it  is  shifted  by  TT . 
Therefore  the  shape  of  the  threshold  structures  is  either 
the  same  as  in  the  radiationless  case  or  it  is  just  turned 
upside  down  as  e.g.  in  RQ  in  fig. 6  at  E,j,  +  W. 

For  E^O.607,  cx»2,N-0,  7 T  and  excitation  of  the  target 

the  quantity  ot(p^.-p^)  is  just  at  the  first  zero  of  the 
Bessel  function  Jg(£(p^-p^))  and  this  explains  the 
accidental  zero  of  Qq  in  fig. 7. 
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4.  Discussions  and  conclusions 
We  have  seen  that  the  low  frequency  approximation  given 
in  (3)  explains  all  qualitative  features  of  resonances  and 
thresholds  in  free-free  transitions  except  the  AC  Stark 
shift. The  quantitative  error  of  the  low  frequency  approxi¬ 
mation  -  corrected  for  the  AC  Stark  shift  -  has  always 
been  below  a  few  percent  (  see  fig. 9  ).But  the  dependence 
of  this  error  on  the  various  parameters  will  probably 
depend  on  the  dimension  of  the  space .Therefore  it  would 
not  be  of  general  interest  to  add  here  a  detailed 
investigation  of  the  numerical  error  of  the  low  frequency 
approximation  in  our  1 -dimensional  model. 


What  advice  can  we  give  to  an  experimentalist , who  wants 
to  utilize  free-free  transitions  in  order  to  look  for 
resonances  and  cusps  in  electron-atom  scattering? 

The  most  interesting  property  of  free-free  transitions 
shown  in  this  paper  is  that  only  rapidly  varying  parts 
of  the  radiationless  scattering  cause!  any  free-free  signal 
in  foreward  direction. This  fact  may  be  used  in  an  experiment 
to  project  out  resonances  and  threshold  effects. The  price 
to  pay  for  this  removal  of  the  background  is  a 
reduction  of  the  number  of  electrons  which  contribute  to 
the  signal. As  shown  in  figs. 2-5  the  ratio  between  the 
resonance  height  in  Tjj,N40  and  in  Trl  is  0.1  in  the  most 


favorable  cases. 


In  an  experiment  it  is  not  possible  to  measure  at  an 
angle  of  exactly  =O.By  proper  adjustment  of  the  laser 
polarization  vector  iT  it  is  possible  to  fulfil  the  relation 

(8) 

.1  ^ 

for  any  choice  of  p^  and  p^,just  turn  £.  perpendicular  to 
the  momentum  transfer  q=pf-p^.At  the  same  time  we  want  to 
see  a  strong  resonance  signal  and  according  to  (5) 
and  g  p^  should  be  big  enough  to  give  a  big  value  for 
the  Bessel  functions. Therefore ,  g  should  be  as  parallel  to 
p^  and  p^  as  possible. All  these  conditions  can  be  met  best 
if  p^  and  p^  are  as  parallel  to  each  other  as  the  construction 
of  the  electron  spectrometer  allows  and  if  g*  is  in  the 
direction  exactly  in  between  and  p^. 

As  indicated  by  the  figs. 2-5  resonance  signals  are 
biggest  and  clearest  in  T+2  or  T_2  at  Ei=ER-  co  or  Ei=ER+ 
respectively  if  the  laser  power  density  is  so  big  that  3?  Pp 
and  <x  pi  are  in  between  1.5  and  2.2  i.e.  if  J  ^Cocp)  is 
near  its  first  maximum. In  this  case  the  resonance  signal 
comes  from  the  term  J^,,(c?pp  f  co  ,  $)  J^(^p^)  in  (5). 

This  means: The  electron  comes  in  shifted  by  one  photons 
energy  away  from  the  resonance  energy. Then  it  absorbs/emits 
one  photon  and  can  go  into  the  resonance  state. The  resonance 
decays  and  the  outgoing  electron  absorbs/emits  another 
photon  on  its  way  out. For  this  choice  of  the  parameters 
the  resonance  signal  in  T^2  is  about  10%  of  the  resonance 
signal  in  the  radiationless  scattering, which  is  the  best  one 
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can  hope  for. In  addition, for  these  values  of  the 
parameters  the  other  resonance  structures  in  T ^  (  i.e. 
the  ones  at  E^,E^+  2U;  ,Ep-t  3^  »etc.)  are  quite  small  and 
this  decreases  the  sources  of  confusion  if  several  elastic 
resonances  occur  close  together. 

For  the  investigation  of  thresholds  choose  a  moderate 

power  density  of  the  laser  so  that  only  first  order 

contributions  give  a  strong  signal  to  the  free-free 

amplitudes  and  look  at  T  ^  or  for  structures  of  the 

length  co  .Two  sharp  drops  in  T+^  at  E&-  cj  and  Eq  and  a 

smooth  behavior  in  between  indicate  a  threshold  at  E  .A 

a 

similar  structure  in  T  „  between  E  and  E  +OJ  indicates 

-1  a  a 

the  same  threshold  at  E  . 

a 

As  indicated  in  figs. 6-8  the  absolute  signal  in  T  N,N$0 
is  extremely  weak  but  the  relative  threshold  effect  is  1 
and  can  therefore  be  seen  clearly. This  is  in  sharp  contrast 
to  elastic  scattering  where  the  relative  threshold  effects 
are  generally  small. Usually  these  effects  are  observed  in 
electron  impact  excitation, where  their  relative  effect  is 
bigger  than  in  the  elastic  channel. Our  idea  is  unusual  in 
so  far  that  it  filters  out  threshold  effects  in  a  process 
which  is  elastic  with  respect  to  the  electron-target 
interaction. The  big  disadvantage  of  our  method  is  the 
extremely  small  absolute  size  of  the  signal  and  there  is 
the  possibility  that  in  a  real  experiment  machine  noise 
will  bu  ry  the  weak  threshold  signal  and  make  it  impossible 
to  utilize  our  idea  in  the  laboratory. 
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In  any  case  for  resonances  or  thresholds  the  photon 
energy  should  be  chosen  larger  than  the  energy  interval 
over  which  the  radiationless  scattering  varies  rapidly  in 
order  to  avoid  confusing  overlaps  of  the  structures  in 
free-free  transitions. 

Can  "accidental"  zeros  like  the  one  in  shown  in  fig. 7 
(  i.e.  those  values  of  parameters  for  which  S(p^-pJ)  is 
at  a  zero  of  some  Bessel  function  )  used  for  anything?We 
don't  think  so, because  their  position  depends  strongly  on 
the  exact  laser  power  and  in  an  experiment  the  laser  power 
varies  in  time  and  space  and  the  electrons , collected  in 
the  detector, have  experienced  quite  different  laser  powers 
during  their  scattering  process .Therefore , such  an  accidental 
zero  would  be  completely  smeared  out  in  any  real  experiment. 

There  is  also  the  interesting  possibility  that  a  free- 
free  experiment  might  be  used  to  probe  the  laser  field 
strength. In  ref. 11  it  has  been  shown  how  the  results  of  a 
measurement  of  jfjj(E^  ffr  )  J  ^  as  function  of  N  can  be  used 
to  calculate  the  average  power  density  of  the  laser  field 
in  which  the  experiment  has  been  performed.lt  is  probably 
possible  to  generalize  the  results  of  ref. 11  and  to  work 
backwards  to  even  more  detailed  information  about  the 


laser  field 
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Appendix 

In  this  section  we  explain  our  model  in  detail  and  show 
how  we  calculate  the  amplitudes  for  all  processes. 

We  work  in  a  1-dimensional  space  and  choose  x  as  space 
coordinate .We  describe  the  electron-target  interaction  by 
a  2-channel  square  well  potential  of  the  form 


v (x)  -  (v-  V22  y  r  -ix'> 

V12=V21  ^  ^  so  tha't  1/  a  self  adjoint  operator. 

@  is  the  unit  step  function. The  matrix  of  excitation  energies  is 


( A2) 


We  cut  the  x-axis  into  the  three  interval  s 
I1  =  (+R,oo)  ,  I2=(-R,+R)  ,  I3=(-oO,-R) 

The  target  is  in  its  ground  state  initially  and  the  electron 

comes  in  with  momentum  k.  .Then  the  matrix  of  the  electron 

in 

momentum  in  intervalls  1^  and  Ij  is 


K-( 


kQ  0 
0  h 


kin 


0 


(A5) 


according  to 


In  I2  we  set  and  construct  the  orthogonal 

matrix  (J/  which  diagonalizes 

QJD  =  £  with 

Vi  0 

0  G*2 

[kin  4  -2b  S]1  is  the  diagonalized  momentum  matrix  in  I2. 


(A4) 


We  assume  a  single  mode  laser  field  and  use  the  dipole 
approximation . 
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The  total  Hamiltonian  H  for  the  motion  of  the  electron 
under  the  simultaneous  influence  of  the  target  and  the 
laser  field  is 

H={-?5-  - £(a++a)  j-  +  (Q)  +  ^(x)  (A5) 

is  an  abbreviation  for  =(  2  7T c^  /o->  L^)3' where  is  the 
quantization  volume  of  the  electromagnetic  field. p  is  the 
electron  momentum  operator, a  and  a+  are  the  annihilation 
and  creation  operators  for  a  laser  photon. We  write  I N  for 
an  eigenstate  of  a+a  with  eigenvalue  N  6- . 

First  let  us  look  at  eigenfunctions  of  H  in  the 
intervalls  I ^  and  1^. There  H  is  diagonal  with  respect  to 
the  target  states. We  neglect  photon  depletion  effects  in 
the  sense  that  we  set  V  N+n"  =  /F  in  the  coupling  strength 
between  various  photon  number  states, where  N  is  the  initial 
number  of  photons  in  the  laser  beam  and  N+n  is  the  number 
of  photons  in  intermediate  or  final  states. Using  the 
recursion  formula  of  the  Bessel  functions 
2n  Jn(y)  -  y^n.qCy)  +  (A6) 

we  see  that  functions  of  the  form 

Cp(kQ,x,N)=  exp  (  ikQx  )  2L  Jn(«  feQ)  |  N+n)>  (A7) 

n 

are  eigenfunctions  of  the  operator 
2  ,  ep 

H_  +  hcja+a  -  —  6  (a+  +  a  )  (A8) 

2m  mc  r 

with  eigenvalue 

Ncu  +  /  2m 


(A9) 


0<  is  an  abbreviation  for  o<=e>(N Jl^/  me  t  u*. 

The  quantity  2  is  the  amplitude  of  the  corresponding 

classical  vector  potential  A  of  the  laser  field. Note  that 
kQ  is  independent  of  n  in  (A7)  because  recoil  effects  of 
the  electron  are  neglected  in  dipole  approximation. In  order 
to  start  from  the  most  general  eigenfunctions  of  H  we 
must  consider  that  the  functions  in  (A?)  are  degenerate 
in  two  ways : 

1. :  (f  remains  an  eigenfunction  of  operator  (A8)  to  the 

same  energy  (A9)  if  we  reverse  kQ  i.e.  if  we  replace 

k0  by  -ko* 

2. : We  get  an  eigenfunction  to  the  same  energy  if  we  replace 


N  by  N+L  and  simultaneously  replace  kQ  by  kL  where 

kf  /  2m  +  Loo  =  k2  /  2m 
L.  o 


Therefore  the  most  general  eigenfunction  of  $={j  with 
eigenvalue  E  in  the  intervall  1^  is  the  2  component 
column  vector 


f 


'  (E,x)  =  Z  [r\Lf(kL>x^+  w1,Lf'kL'x'LA 
1  *  L  \r2,L  fhL»x^+  w2,Lf(‘hL’x’LHy 


(A10) 


where  all  kL  and  h^  are  given  by 
E=  k2  /  2m+(tf-H)aJ  =  h|-/  2m+(N+I^J+ 


(All) 


r^  L  and  w^  L  are  arbitrary  complex  constants  to  be  fixed 
later  by  boundary  conditions. 

In  a  completely  analoguous  way  we  find  for  the  most 
general  eigenfunction  of  ^with  eigenvalue  E  in  interval  1^ 


S  (E,x)-  I  (  Jl.L  <f’(kL'x'LtN>  *  *1,1,  <f<-kL’x'L*N»\ 

-*5  ’  L  V  J2.L  <f^*L,x'I,+^  +  *2,1  ft-hL>x-W 

ji  ^  and  t^  ^  are  again  constants  to  be  fixed  later  by 
boundary  conditions. 

In  the  interval  Ig  (f“^is  not  diagonal  with  respect 
to  t)8P^Gt<  8t/8t/6S  tut  i y*MU-  W  is  and  therefore  we 

'S/ 

first  construct  eigenfunctions  to  H  .In  the  same  way  as 

<v 

before  we  find  for  the  most  general  eigenfunction  of 
with  eigenvalue  E 

X  (EiX)a  £  f8l'L  lf(AL’x'L+M,+  b1,L?C'^'x»L+f,)\ 

*  L  V  a2,L  tf(/Vx’L+^+  b2,Lf  (-/L»x’L+7 

where  /  2m+^(  +  L)to  +  *E«  /  2m*-(N+L)c*>+ 

with  and  given  in  (A4). 

a^  L  and  b^  L  are  again  constants  to  be  fixed  later  by 

boundary  conditions.  ^ 

Because  of  Hy%  -y>H%  -  e  u$2 

/■>/ 

we  see  that  31  ^$2  an  e^6enf unction  of  (HI  with 


(A12) 


(A13) 

(A14) 


eigenvalue  E  in  interval  I2. 

At  x*+oo  we  choose  the  boundary  condition  that  there  is 
only  one  incoming  wave  with  the  target  in  state  1  and  the 
laser  and  electron  in  state  <p(-k^n,x,N)  i.e. 

wi,L  "  ^i,1  '  ^L,0  (A15) 

At  x—oo  we  choose  the  boundary  condition  that  there  is 

no  incoming  wave  at  all  from  the  left  i.e. 

T*0  for  all  i  and  all  L  (A16) 

At  x»+R  and  x»-R  we  require  that  the  wavefunction  ^  is 
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continuous  and  has  a  continuous  first  derivative. This 
gives  A  equations  for  2-component  column  functions  or  8 
equations  for  linear  combinations  of  photon  number  states. 
These  equations  must  be  fulfilled  for  the  coefficients 
of  each  photon  state  separately  and  therefore  each  photon 
state  gives  us  8  equations  which  connect  the  constants 
ai  L,bi  L’^i  L,ri  L  ^  i=>1  ).Some  of  these  equations 
are  inhomogeneous. 

Al  together, we  get  an  infinite  inhomogeneous  system  of 
linear  equations  to  determine  the  free  constants. In  our 
model  calculations  we  could  only  handle  a  finite  number 
of  equations  with  a  finite  number  of  unknowns  and  therefore 
we  had  to  cut  off  the  system  of  equations. Out  of  the 
coefficient  matrix  of  the  system  of  linear  equations  we 
have  cut  out  the  (21+1 ) • 8  * (21+1 ) • 8  matrix  centered  at 
the  8x8  block  which  comes  from  the  L=0  terms  in  (A10,A12, 

A13)  and  the  n=0  terms  in  (A7) .Accordingly  we  have  cut 
out  of  the  inhomogenity  vector  a  (21+1)8  component  piece 
centered  at  the  block  which  comes  from  the  n=0  terms  in 
(A7).Then  we  have  solved  these  (21+1)8  coupled  linear 
equations. 

We  have  found  rapid  convergence  of  the  results  with 
increasing  1  values  as  soon  as  l>|&p[.For  all  calculations 
shown  in  figs. 2-9  it  has  been  sufficient  to  choose  a  1 
value  between  5  and  10. The  rapid  convergence  can  be  understood 
from  the  fact  that  Bessel  functions  decrease  rapidly  as 
soon  as  the  absolute  value  of  the  order  becomes  larger  than 


the  absolute  value  of  the  argument. If  we  take  a  larger 
value  of  l,then  we  include  coefficients  which  contain 
Bessel  functions  of  higher  orders. 


AS  a  last  step  we  calculate  the  quantities 
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These  are  the  quantities  plotted  in  figs. 2-9  as  function 
of  the  incoming  electron  energy 

All  these  calculations  can  also  be  performed  for  the 
case  that  the  target  is  in  state  2  initially .The  only 
change  is  to  replace  OC  of  (A3)  by 
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DC 
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and  to  replace  (A15)  by 


wi,L 


i,2  ®L,0 
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Figure  captions 

Fig.1 :Plot  of  the  diagonal  elements  of  the  model  potential 
for  electron-target  interaction  used  in  all  numerical 
calculations. The  offdiagonal  elements  vi2=V21  ^  not  Plott:ed 
in  fig.1  )  are  of  the  same  shape  as  but  of  different 
depth. Bound  states  of  the  1-channel  potential  V22+E12 
become  Feshbach  type  resonances  of  the  full  coupled 
2-channel  potential. Their  width  P  is  determined  by  the 
magnitude  of  V^^or  small  we  find  T  ©C  (  )2.See 

the  two  resonances  in  the  upper  lines  of  fig. 2  and  fig. 4. 

A  threshold  effect  occurs  at  the  rim  of  the  upper  well 
i.e.  at  £=£^2* 

Fig.2:Plots  of  resonance  structures  in  radiationless 
scattering  and  in  free-free  transit ions. For  more 
explanations  see  main  text. 

Fjg.3:Plot8  of  resonance  structures  in  free-free 
transitions. For  more  explanations  see  main  text. 

Fig.4:Plot8  of  resonance  structures  in  radiationless 
scattering  and  in  free-free  transitions. For  more 
explanations  see  main  text. 

Fig.5:Plots  of  resonance  structures  in  free-free 
transitions. For  more  explanations  see  main  text. 

Fig.6:Plots  of  threshold  structures  in  radiationless 
scattering  and  in  free-free  transitions. For  more 
explanations  see  main  text. 

Fig.7sPlot8  of  threshold  structures  in  free-free 
transitions. For  more  explanations  see  main  text. 
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Fig.8:Plots  of  threshold  structures  in  free-free 
transitions. For  more  explanations  see  main  text. 

Fig.9:Some  examples  for  the  AC  Stark  shift  of  resonances 
in  free-free  transit ions. The  solid  lines  are  the  results 
of  numerical  calculations. The  broken  lines  are  the  results 
of  eq.(3)»For  more  explanations  see  main  text. 
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